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Abstract 



In the framework of chiral perturbation theory we investigate whether, in the presence of 
nucleons, a many-particle system may lower its energy density by changing its mesonic vac- 
uum from a homogeneous to a spiral configuration. Numerical values for the critical particle 
number densities for the spiral configuration to occur are given. They lie in the region of nu- 
clear matter density where the applicability of chiral perturbation theory should be checked 
carefully. 
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Introduction 



Quantum chromodynamics, the field theory of the strong interaction, is an asymptotically 
free theory, i.e. the coupling constant decreases as the energy increases. On the other hand, 
at low energies the coupling constant becomes too large to make an expansion in it; QCD 
becomes non-perturbative. Besides lattice calculations, the effective field theory is the only 
systematic approach to deal with the low-energy dynamics of QCD. The symmetry of the 
QCD Lagrangian in the chiral limit (i.e. at vanishing quark masses) is supposed to be broken 
spontaneously leading to massless Goldstone bosons by virtue of the Goldstone theorem. 
Depending on whether the s-quark is considered as light or not, these Goldstone bosons may 
be identified either as the pseudoscalar octet or as the three pions. Since chiral symmetry is 
explicitly broken by non-zero quark masses the Goldstone bosons are not exactly massless but 
anyhow relatively light compared to 1 GeV, the typical scale of QCD. At low energies they are 
the fundamental degrees of freedom, whose dynamics is dictated by symmetry constraints. 
This leads to the effective field theory of QCD, chiral perturbation theory, formulated in its 
modern form by Weinberg, Gasser and Leutwyler. 

Here, we consider the s-quark as heavy and deal with pions only. In the pure mesonic 
sector the ground state is a constant pionic field configuration. As we will discuss, the pions 
may be represented by a four-dimensional unit vector, such that the ground state just stands 
for a constant vector field. However, chiral perturbation theory may be extended to include 
also baryons (nucleons in the case where the s-quark is considered as heavy). This makes 
it possible to effectively describe the interactions between one nucleon and pions. It turns 
out that the fermion may lower its energy if the vector field representing the pions is not 
constant but varies throughout space. This behavior increases the pionic energy and for low 
nucleon densities this rise of the pionic energy outweighs the decrease of the fermionic energy. 
However, if we do not consider one nucleon only but a many-particle system (neglecting 
interactions between the fermions) and if the particle number density of nucleons exceeds 
some critical value the dropping of the fermionic energy becomes so large that the many- 
particle system may lower its total energy density by changing the pionic configuration from 
the homogeneous to an inhomogeneous phase. 

In this thesis we investigate whether it is possible that the vector field representing the 
pions could show a spiral phase in the ground state if we take nucleons into account. The 
physical parameter we can control is the density of nucleons. In the framework of chiral 
perturbation theory we calculate the energy of a spiral configuration (to be introduced in 
section I2.2p and try to find out what critical value the fermionic particle number density 
must exceed for the spiral phase instead of the homogeneous phase to be favored. The 
spiral configuration we consider is such that the nucleons see a constant background^, which 
simplifies the solution of the Dirac equation. 

The idea of this thesis actually originates from the theoretical description of quantum 
antiferromagnets. Antiferromagnets doped with either electrons or holes are nowadays the 

1 Technically speaking they couple to a constant connection and vielbein, see 11.8.21 
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only known approach to high temperature superconductors and have therefore been studied 
also from the theoretical point of view intensely. There are different models for the description 
of quantum antiferromagnets on the microscopic level. The most important are the so called 
Hubbard and t — J model. However, currently these models can be simulated numerically 
only if they contain at most one single hole or electron. Additionally, due to severe fermion 
sign problems and other complications it is not possible to find the exact ground states of 
the two doped models, neither analytically nor numerically. 

For these reasons one constructs a low energy (low temperature) effective theory for 
antiferromagnets, i.e. as in the case of QCD one considers only the low energy regime and 
the degrees of freedom appearing there. At low energies, the spin symmetry of the microscopic 
model is spontaneously broken giving rise to two massless Goldstone bosons (magnons) that 
can be represented by a three dimensional unit vector e(x), the so called local staggered 
magnetization. An effective theory in terms of e(x) for the pure antiferromagnet, i.e. for the 
magnon sector without holes or electrons, was developed by Hasenfratz and Niedermayer in 
|34j . Considering the effective action therein it becomes clear that the ground state of the 
antiferromagnet at low energies corresponds to a constant staggered magnetization. 

Like chiral perturbation theory may be extended to include baryons, doped holes (|35j) 
and doped electrons ( [36J ) in the antiferromagnet may also be described effectively. In [37] 
it was shown that through an inhomogeneous staggered magnetization, namely a spiral con- 
figuration, (which costs energy in the magnon sector) doped holes can lower their energy in 
such a way that the total energy of the system may get lower than in the case of a constant 
staggered magnetization. Depending on the values of the low energy constants and the den- 
sity of holes, spiral configurations are hence a certain category of configurations that may 
then describe the ground state. 

Due to the similarities of the two effective field theories, one for QCD and one for anti- 
ferromagnets, it was natural to ask whether the corresponding phenomenon to spiral phases 
observed in antiferromagnets could occur also in QCD. 

The answer to this question is positive. However, as we will show, the critical particle 
number density for the spiral phase to occur will turn out to lie in the region of nuclear 
matter density. It is anything but obvious whether chiral perturbation theory remains valid 
at such high densities. In order to check the reliability of our results we investigate higher 
orders (in different ways) and calculate the corrections coming from these. If they are not 
too large we might get a reason to trust the results and to believe that the qualitative picture 
may be correct. Due to these uncertainties we will not calculate errors for the critical number 
densities originating from errors of the low energy constants. The calculations are performed 
using the low energy constants obtained in the chiral limit, but in order to analyze the 
sensitivity of the results with respect to higher orders we also do the calculations with the 
physical low energy constants. 

The presence of a spiral configuration could be viewed as the phenomenon of pion conden- 
sation. This is supposed to occur in superdense matter like e.g. in nuclear matter or inside a 
neutron star. Spiral configurations in the mesonic vacuum could therefore be an alternative 
approach towards such phenomena. 



Chapter 1 



Some Remarks on Chiral 
Perturbation Theory 

This chapter is intended to give a brief summary of chiral perturbation theory as the effective 
theory of QCD. More detailed introductions can be found in [U [5j El [12j fl5| f!6] . 

1.1 The Need for an Effective Field Theory 

The Lagrangian of quantum chromodynamics (QCD) reads 

6 1 
£qcd = £ 9f{0 ~ m f )q f - -^Q^aQt- (1-1) 

/=i 

Here / is the flavor index, qf a triplet of spinors (red, green and blue quarks), 

8 A 

a=l 

is the covariant derivative with the eight gluon fields Au,aj the strong coupling constant g 
and the Gell-Mann matrices A a , and 

Gfj,i/,a — dfj,A-u,a du-Afi,a 9 fabc*A[j,,b-Av,c (1"^) 

is the field tensor for the vector potentials, i.e. the kinetic term for the gluons with the SU (3) 
structure constants f a b c - We have adopted Feynman's notation Jj) = 7 M -D M . 

Since QCD is a non-Abelian theory (f a bc 7^ 0) the gluons carry the same color charge as 
the quarks and thus may interact with each other. The structure of QCD is therefore much 
more complicated than that of QED and, what is more important, due to the non-Abelian 
nature QCD is an asymptotically free theory, i.e. the (3 function /3 = fidg/dfi (where fj, is the 
renormalization scale, i.e. the typical energy) is negative^] This means that, in contrast to 
QED, the coupling constant decreases with increasing energy. Since at low energies the strong 
coupling constant becomes quite large, perturbation theory, i.e. an expansion in the coupling 
constant, is possible only for high energies but will fail in the low energy regime. Only two 
rigorous approaches are known to treat the strong interaction at low energies: Lattice QCD 
and the effective field theory called chiral perturbation theory (ChPT). 

1 One can show (see e.g. [SS] or [35]) that the fi function of QCD would be positive only if the number of 
quark flavors exceeded sixteen. 
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1.1. The Need for an Effective Field Theory 



1.1.1 Effective Field Theories 

Effective field theories (EFTs) are a very important tool in many areas of physics. The 
motivation for an EFT is that one does not have to know a theory at all energy scales in 
order to precisely describe a physical system. In general, EFTs are low energy approximations 
of a fundamental theory and allow an expansion of physical quantities in terms of p/A where 
p denotes momenta or masses being smaller than some scale A (therefore an expansion is 
possible). The basis of EFTs is a statement given by Weinberg and proved by Leutwyler 
([9]) as well as d'Hoker and Weinberg ([E]), which says that a quantum field theory has no 
content beyond unitarity, analyticity, cluster decomposition and symmetries. 

Unitarity assures that X^/ K/I^l*)! 2 = 1 (/• fi na l state, i: initial state) and analyticity 
guarantees causality. The cluster decomposition theorem guarantees locality and states, in 
few words, that two experiments done in sufficiently separated regions of space-time do not 
correlate. Finally, by symmetry we mean Poincare invariance, C, P, T symmetries and 
internal symmetries like e.g. isospin symmetry. 

Weinberg's statement means that a perturbative description in terms of the most general 
effective Lagrangian containing all possible terms compatible with assumed symmetry prin- 
ciples yields the most general S matrix consistent with the fundamental principles mentioned 
above. Hence, an effective Lagrangian is constructed by writing down all possible terms al- 
lowed by the symmetry properties of the underlying theory. Of course, this will be an infinite 
number of terms and we need some method of assessing the importance of diagrams generated 
by the interaction terms of this Lagrangian. In ChPT this will be Weinberg's power counting 
scheme (see section 11.5. lj) . 

There are different types of EFTs (see e.g. [Tj) depending on the structure of the transition 
from the fundamental (energies > A) to the effective (energies < A) level: 

• Decoupling: If we have a theory in which there are light and heavy degrees of freedom 
which are well separated by a scale A and if we only consider energies well below 
A, then the heavy degrees of freedom will not matter. We can integrate them out 
of the generating functional and get a Lagrangian with the light degrees of freedom 
only (Appelquist-Carazzone theorem, see [22 J). This effective Lagrangian will contain a 
renormalizable part and non-renormalizable couplings which are suppressed by inverse 
powers of A. An example for such an EFT might be the Standard Model (as an EFT 
of a yet unknown fundamental theory). 

• Non-decoupling: Due to some phase transition the degrees of freedom at low energy 
differ from the ones of the underlying theory. This is the case in QCD where the phase 
transition is spontaneous symmetry breaking of chiral symmetry. In the low energy 
regime we only deal with color-neutral hadronic states which will be the Goldstone 
bosons in the mesonic sector (originating from spontaneous symmetry breaking) and 
nucleons in the baryonic sector. An effective Lagrangian of that sort will not be renor- 
malizable in the classical sense. The reason is that it contains infinitely many terms and 
due to the nonlinear transformation properties of the Goldstone bosons even the lowest 
order effective Lagrangian contains infinitely many fields. To get a dimensionless action 
we thus need coupling constants of negative mass dimensionality. But such interactions 
are not renormalizable. However, this is not regarded as a serious problem since the 
infinities arising from loops in a given order of momentum expansion are absorbed in a 
renormalization of (a finite number of) coefficients in the Lagrangian. So the effective 
Lagrangian is renormalizable order by order. 

Going from the fundamental to an effective level we lose track of the correct high energy 
behavior of the theory. An EFT will yield wrong results if the energy is too high and e.g. 
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heavier degrees of freedom are produced. Thus the domain of utility of an EFT is necessarily 
bounded from above in the energy scale. 

The effective coupling constants, i.e. the coefficients appearing in the effective Lagrangian, 
cannot be determined by symmetry considerations alone but should in principle be calculated 
from the fundamental theory. In the case of ChPT we cannot (yet) solve the fundamental 
theory (QCD) and the effective couplings stay as free parameters which have to be determined 
by some other means, e.g. matched by experiments. 

1.2 Symmetries of the QCD Lagrangian 

Since symmetry and symmetry breaking are the basic concepts in the construction of an 
effective Lagrangian we shall study the symmetries of the QCD Lagrangian in some detail. 

1.2.1 The Chiral Limit 

Due to the confinement of quarks into hadronic states the notion of quark masses is rather 
complicated and needs some extended discussion. However, without going into details, e.g. 
the Particle Physics Booklet [24J tells us that there are three light (u,d and s) and three 
heavy (c,b and t) quarks. The masses of the light quarks are quite below 1 GeV and the 
heavy quarks exceed this scale. 1 GeV is the typical scale of QCD and is associated with the 
masses of the lightest hadrons containing light quarks that are not Goldstone bosons (e.g. 
m p = 770 MeV). For spontaneous symmetry breaking the typical scale is 47ri ? 7r ~ 1.2 GeV 
(see e.g. [21])) which is of the same order. If we constrain ourselves to energies below 1 GeV 
there will never appear any states containing heavy quarks and we can ignore them in the 
Lagrangian, i.e. we set their masses to infinity. 

If we compare e.g. the mass of a proton (m p = 938 MeV) with the sum of two up and 
one down quark we see that the quark masses contribute almost nothing and that the proton 
gets its mass through another complicated mechanism. This is the motivation to send the 
light quark masses to zero as a first approximation. This is the so-called chiral limit and the 
Lagrangian (jl.ip then reads 

£°qcd= E «fifaf-\s m a9?'. (1-4) 

f=u,d,s 

Spinors of massless or ultrarelativistic fermions are eigenstates of the chirality operator 
75 = ^7°7 1 7 2 7 3 , where particles with eigenvalue +1 are called right-handed and those with 
eigenvalue —1 are called left-handed. In this extreme relativistic case chirality equals to helic- 
ity. Defining the projection operators Pr = ^(1 + 75) and Pl = |(1 — 75) we can decompose 
an arbitrary spinor into its right- and left-handed parts which then are eigenstates of 75. 
Thus, with qn t L = Pr,lQ and q Ry L = qPl,r we can then write (jl.4p as 

£qcd= E {W^RJ + QL,f0qL,f)-\Q^,aQ^- (1.5) 

f=u,d,s 

We see that the right- and left-handed quark fields have completely decoupled. This is not 
the case if we have a mass term which mixes right- and left-handed fields. Besides Lorentz 
invariance, SU(3) C gauge invariance and P, T and C symmetries, the Lagrangian fjl .5j) shows 
also a chiral symmetry: It is invariant under global U(3)l x U(3)r flavor transformations 
(QL, Qr) -> (LqL, Rqr) and (q~L,q~R) -> {qLLt , q R P3) where L,R£ U(3) l ,r and q L>R are flavor 
triplets. Since we can make the decomposition U(N) = SU(N) x U(l), we see that £q CD 
shows a global SU(3)l x SU(3)r x U(1)l x U(1)r symmetry. 
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1.2.2 Symmetry Currents 

Prom these symmetries we get several conserved currents due to Noether's theorem (see e.g. 
[5]). The U (1)_r,l symmetries will lead to the currents <Jr > l'Y^QR,L' By adding and subtracting 
respectively these currents we get a singlet vector and a singlet axialvector current (according 
to to their behavior under parity transformation): 

V» = q R ^q R + qil^QL = q^q (1.6) 
A" = qR^qR - qLl^qi = qi^im (1.7) 

The current can be interpreted as total current when doing the U(1)l and U(1)r trans- 
formation with the same phase, while comes from transformations with opposite phases 
(if a U(l) transformation with phase 6 yields the current J^, then a transformation with 
phase kQ (k £ R) leads to kJ^). 

From the conserved charge is given by Qy = J d 3 xV° = f d 3 xq^q which counts 
the number of quarks minus antiquaries (see e.g. [6]). Thus Qv/3 is the number B of 
baryons and the U(l)y symmetry simply expresses baryon number conservation and leads to 
a classification of hadrons into mesons (B = 0) and baryons (B = 1). 

The current is conserved only at the classical level. This symmetry is broken by 
quantization due to anomalies (see e.g. [21j). 

The interesting symmetry is SU(3)l x SU(3)r which leads to eight left-handed and eight 
right-handed conserved currents: 

L^ a = q L ^Y qL ^ 
R^ a = qia^QR- (1-9) 
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The charge operators read 



Q a L = j d 3 xL°> a = J d 3 xq[^-q L (1.10) 
Q% = j d 3 xR°> a = f d 3 xq j R ^-q R (1.11) 



and are the generators of SU(3)l x SU(3)r since they satisfy the commutation relations 
corresponding to the Lie algebra of this group. Again one considers the vector and axialvector 
linear combinations 

y^a = R ^a + L ^a = q (1.12) 

A »,a = R H,a_ £/^ = g 7 M 75 *L Qi (1.13) 

as well as 



Qv = Qr + Ql (1-14) 

Q a A = Q a R-Qh (1-15) 

Note that the charge operators Q v satisfy the commutation relations for SU(3)y, whereas 
Q a A do not form a closed algebra, i.e. the commutator of two axial charge operators is not 
again an axial charge operator and there is nothing like SU{3)a- 
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1.3 Spontaneous Symmetry Breaking 

Spontaneous symmetry breaking takes place if the Lagrangian of a system has a given sym- 
metry but the ground state is not invariant under that symmetry; rather the system has 
chosen one of (possibly many) ground states, all related by symmetry transformations. A 
classical example is a ferromagnet described by spin-spin interaction, whose Hamiltonian is 
invariant under rotations. However, below the critical temperature the ground state is one 
in which all the spins point in the same direction, which is surely not rotationally invariant. 
The direction of the spins is random (in the absence of an external magnetic field), hence the 
system has infinitely many degenerated ground states. 

As we saw that £q CD is invariant under SU(3)l x SU(3)r x U{l)y, it is interesting to 
ask whether the ground state of QCD is invariant under the same symmetry group. To study 
this, let be an eigenstate of Hqqd with some energy E and with positive parity, i.e. 
H^ CD \ij) = E\ip) and P\ip) = Then the state Q A \ip) is also an eigenstate of Hqqd 

with the same eigenvalue E but with opposite parity: 

#qcd w> = q\h%cM = norm 

PQ a A \i>) = PQ a P- 1 P\iP) = -QHV), 

■ since Q a A is time independent and PQ a A P~ l = —Q a A 



where we have used 



^QCD> Q°A 
since Q A is an axialvector. 

Consequently, for any state of positive parity one would expect the existence of a degen- 
erate state of negative parity. However, this is not observed in the low-energy spectrum of 
hadrons. Though, for these arguments we (not obviously) assumed that the axial charges 
Q a A annihilate the ground state, i.e. Q a A \Q) = (see [5]). But since there is no experimental 
evidence for parity doubling one concludes that the Q a A do not annihilate the vacuum, i.e. 
the ground state of QCD is not invariant under SU(3)l x SU(3)r x U(l)y. In [23] it was 
shown that in the chiral limit the ground state must be invariant under SU(3)y x U(l)y, 
i.e. the vector charges Q v and Qy annihilate the ground state; <3y|0) = Qy|0) = 0. Thus, 
while £q CD is invariant under SU(3)l x SU(3)r x U(l)y, the ground state is invariant under 
SU(3)y x U(l)y only, which is spontaneous symmetry breaking. 

1.3.1 Goldstone's Theorem 

Goldstone theorem states that for every spontaneously broken continuous symmetry, i.e. for 
every generator that does not annihilate the ground state, there is one particle of zero mass 
and spin in the spectrum. These are called Goldstone bosons. 

Let us consider a theory involving several fields (j)i(x), described by a Lagrangian of the 
form C = (derivatives) — V(4>). Suppose that the potential V(<j>) has a set of degenerate 
ground states and that the system has chosen spontaneously one of them, let us call it 0o- 
Expanding V about this minimum yields 

V{$) = V{^) + \{$-^)S-^) j ^- +..., (1.16) 
where the linear terms disappear since 4>q is a minimum. The coefficients 



d 2 v 



d<fiid(f)j 



(f>=tj>0 



m% (1.17) 



define a symmetric matrix whose eigenvalues give the squared masses of the particle excita- 
tions above the ground state (vacuum) (pQ. This can be seen by performing an orthogonal 
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transformation of the fields such that the mass matrix gets diagonal. Since mf - is a positive 

semidefinite matrix (as (fio is a minimum) these eigenvalues cannot be negative. We show 
now, that for every continuous symmetry of C which is not a symmetry of 4>q the matrix 
mf - yields a zero eigenvalue and we therefore have a massless particle. Suppose that C is 

invariant under a symmetry group G of order uq (thus G having uq generators) while (j)Q is 
invariant under a subgroup H C G of order nu only. Under G an infinitesimal symmetry 
transformation of the fields takes the form 



% = <t>i 



+ ie a T^j, (1.18) 



where T a are the representation matrices of the generators of G. Since the potential V{4>) is 
invariant under G, we have 

- dV 
V(<f>) = V{4>') = V(ct> + ie a T a 4>) = V{fa + ieaTftfafa + ie a T^j, ...) = V(<f>) + —%t a T^ 

0(pi 

which yields 

dV 

g^ieaT^j = 0. (1.19) 
Differentiating with respect to 4>k and evaluating the resulting expression at = <j)Q we get 



d 2 v 



%oq> k 



4>=4>o 



dV 

iealZfaj + 7TT r r = =► m^T^oj = 0, (1.20) 



4>=<f>o 




since the equation holds for every € a . This can be written more compactly as Tn^(T a (pQ^i — 0. 
This equation is trivially fulfilled if a £ {1, 2, }, i.e. if T a is a generator of the subgroup 
H, because 4>q is invariant under H and hence T a (f)Q = 0. But for a £ {nn + 1, •••,«g}> i- e - 
for the generators belonging to G but not to H, we have T a (pQ ^ 0, since cfto is not invariant 
under the full group G. Then the above equation tells us that T a 0o is an eigenvector of mj^ 
with eigenvalue zero. Thus the matrix has uq — fin zero eigenvalues and therefore there are 
tig — n H massless particles (Goldstone bosons) in the spectrum. These Goldstone bosons 
carry the same quantum numbers as the generators that do not annihilate the vacuum. 
Furthermore, it follows (see [5]) that the matrix element of the symmetry currents J^(x) that 
lead to the generators of G not annihilating the ground state between the vacuum and the 
massless one-particle states |-7r a ) is non-zero, 

(0|^(0)|vr a )/0. (1.21) 

Using this, it can be shown that, in the chiral limit, the Goldstone bosons do not interact 
with each other at zero momentum (see. e.g. [16] or [7]). This feature is essential for the 
consistency of ChPT; otherwise the power counting of ChPT (see 11.5. ip would break down. 



1.3.2 The Scalar Quark Condensate 

Remembering the example of the ferromagnet below the critical temperature, we see that, as 
soon as the symmetry of the system has broken down spontaneously, the magnetization (M) 
becomes non-zero. (M) can therefore be identified as an order parameter of the system. In a 
similar maimer we can detect a spontaneous breakdown of chiral symmetry by investigating 
the scalar quark condensate (qq) = (0\qq\0) which is the order parameter of QCD (see e.g. 
|7J). A non-zero value of this order parameter in the chiral limit is a sufficient (but not 
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necessary) condition for spontaneous symmetry breaking in QCD (see [5] or [H]). If the 
axial charges Q a A do not annihilate the vacuum it is consistent to assume that (qq) ^ (see 

e.g. mi). 

As we have 8 generators Q a A which do not annihilate the ground state, Goldstone's theorem 
now allows us to conclude that there must be 8 massless pseudoscalar particles \n a ) whose 
coupling to the currents is non-vanishing: 

(0\A«(x)\7r\p)) = ip ,F 5 ab e-^. (1.22) 

This form of the matrix element comes from symmetry considerations: Due to translational 
invariance we have (0|^4^(x)|7r b ) = e~ tpx (0\ A^(0)|7r 6 ) and the only quantity in question bear- 
ing a Lorenz index is the momentum p^ of the pseudoscalar state. The constant of propor- 
tionality Fq is called the pion-decay constant in the chiral limit. It measures the strength 
with which the Goldstone boson \n b ) decays via the axial vector current into the hadronic 
vacuum. 

Acting in (jl.22p on both sides with <9 M we conclude that, in the chiral limit (where d^A^ = 
0), the Goldstone bosons are indeed massless: 

= (0\d^A a ^(x)\TT b ( P )} = p 2 F 5 ab e- ipx p 2 = m 2 = 0. (1.23) 

I. 3.3 Goldstone Bosons of QCD 

Since in the chiral limit of QCD we suppose a spontaneous symmetry breakdown from 
SU(3) L x SU(3) R x U(l) v to SU(3) V x U(l) v and since SU(3) has eight generators, we 
expect 8 massless Goldstone bosons. Their symmetry properties are tightly connected to the 
generators which are responsible for them, i.e. to those which do not annihilate the vacuum. 
In QCD this are, as seen above, the axial generators Q A . Thus the Goldstone bosons are, 
as already mentioned, expected to transform with a negative sign under parity and hence to 
be pseudoscalar s. Indeed, the spectrum of QCD shows eight candidates for the Goldstone 
bosons, namely the pseudoscalar octet (tt,K,ti) = (ir°, tt^, K°, K°, K^, rj). These particles 
are of course not massless but this is interpreted as a consequence of the explicit breaking of 
chiral symmetry due to finite u-, d- and s-quark masses in the Lagrangian (jl.ip . E.g. in [5] 
the resulting masses of the Goldstone bosons are calculated (the results are quoted in section 

II. 4.2p . Nonetheless, the Goldstone bosons of QCD are much lighter than all other hadrons. 
This statement becomes even more accurate if we restrict ourselves only to u- and d-quarks 
and neglect the s-quark (which is considerably heavier than the u- and d-quarks) in the La- 
grangian, i.e. we only consider energies where states including s-quarks never appear. Then 
£q CD shows an SU (2)l x SU (2)r x U (l)y symmetry which breaks down to SU (2)y x U (l)y, 
resulting in three Goldstone bosons, as SU{2) has three generators. These are then the pions 
which, since the u- and d-quarks are indeed very light and thus the chiral limit is a better 
approximation than in the case of SU{3), show strikingly small masses (about 139 MeV for 
vr ± and about 135 MeV for vr°). 

1.4 Construction of the Effective Lagrangian 

Since in QCD at low energies the Goldstone bosons are the only degrees of freedom, this are 
the objects to be contained in an effective Lagrangian. 

1.4.1 Representation of the Goldstone Bosons 

Let us denote the symmetry group of the Lagrangian by G and the smaller symmetry group 
of the ground state by H. Due to spontaneous symmetry breaking we have n = uq — njj 
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1.4. Construction of the Effective Lagrangian 



Goldstone bosons fa which we collect in the vector <3? = ifa, ...,cp n ). The set of all these 
vectors is a vector space, let us call it M. We can define the action of the symmetry group 
G on M by a mapping / : G x M — > M with 

$^$' = fig,$), geG. (1.24) 

This defines an operation of G on M if the mapping / satisfies 

/(e,$) = l V$ (1.25) 
f( gi ,f(g 2 ,$)) = f( gi g 2 ,$) V 9 i, ?2 €G,V$, (1.26) 

where e is the identity of G. Since the subgroup H leaves the ground state <1> = invariant 
we have f(h, 0) = for all h G H and hence 

f(gh, 0) = f{g, f(h, 0)) = f{g, 0) \/g G G, h G H. (1.27) 

This shows that for all elements of a given coset gH = {gh\h G H} oi g, f maps the origin 
onto the same vector in M. Furthermore, this mapping is injective; one simply shows that 
for two elements g and g' of G with g' 4 gH it follows /(<?, 0) ^ f(g',0) (see e.g. [5])- So the 
mapping is bijective on the image of f(g, 0) (as it is automatically surjective on the image) 
and can there be inverted. Hence, to each <3? corresponds a coset gH with appropriate g 
such that <3? = f[gh,0). We have therefore found an isomorphic mapping between the set 
{gH\g G G} of all cosets and the Goldstone boson fields. The set of all cosets is called 
quotient G/H and therefore the Goldstone boson fields are said to live on the coset space 
G/H. 

The transformation behavior of a $ = f(gh,0) under a group element g € G is easily 
found by applying f(g,-): 

$' = fig, $) = fig, figh, 0)) = f(g~gh, 0) (1.28) 

The coset corresponding to <&' is thus obtained by multiplication of the coset corresponding 
to $ with g. 

Now, in the case of QCD we have G = SU(N) x SU(N) = {(L,R)\L,R G SU(N)} 
and H = {(V,V)\V G ST/(iV)} which is isomorphic to SU(N) V , with AT = 3 or JV = 2 
depending on whether we include the s-quark or not. The coset of g = (L, R) G G is 
gH = {iLV, RV)\V G SU(N)} and since we have 

(LV, RV) = (LV,RDiV) = (l,RL ] ) (LV,LV), 

V v ' 

it can be written as gH = (1, RL')H. Hence the coset and therefore the Goldstone bosons 
may be uniquely characterized through the SUiN) matrix U = RU. Its transformation 
behavior under an element g = (L, R) G G is, as seen above, obtained by multiplication of 
the coset with g: 

g gH = (L, RRlJ)H = (1, RRL ] L ] ) (L, L) H = (1, R iRD) L ] )H, 
=$ eH u 

which yields 

U RUL ] . (1.29) 
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This mapping defines a nonlinear realization of SU(N)l x SU(N)r. One possible and con- 
venient representation of the matrix U is the exponential representation: 

U(x) = exp (i^) = exp (i^) , (1.30) 

where t a are the generators of SU (N) and Fq will turn out to be the pion decay constant in 
the chiral limit (see section ["l,4.2p . In terms of the physical fields we have for SU(3) 



<t>a A a 



/tt + 4=tj \/2vr+ V^+A 
V2tt- -tt° + 4=r? v^i^ 



(1.31) 



V5V 



and for 517(2) 



vr° v / 2vr+ 



-TT 



1.4.2 The Effective Lagrangian 

As mentioned in section 11.1. 1| the effective Lagrangian contains all possible terms allowed 
by the symmetry properties of the underlying theory. For our case this means that the 
Lagrangian should be invariant under SU(N)l x SU(N)r x f7(l)y with N = 2, 3, while the 
ground state should only be invariant under SU(N)y X U(l)y. Incorporating the SU(N) 
matrix U it contains eight (three) pseudoscalar degrees of freedom, which, under SU(N)y 
transform as an octet (triplet). This corresponds to the observations in nature as the hadrons 
organize themselves into approximately degenerate multiplets transforming under irreducible 
representations of SU(N). 

It is easily checked that the ground state (j){x) = and thus Uq = 1 is invariant under 
SU (N)v, since in this case we have R = L = V and (|1.29p yields 

U' Q = VU V ] = VV ] = 1 = U Q . 

On the other hand Uq is not invariant under axial transformations, i.e. for L = A and R = A' 
(corresponding to the case of opposite phases, see section fl.2.2j) : 

U' Q = A^U Q A^ = A^A^ ^t = U , 

just as we expect it due to spontaneous symmetry breaking. The transformation behavior 
of the Goldstone boson fields cp a may be verified to be that of a multiplet under SU(N)y. 
Thus, we can start writing down all possible allowed terms containing U, W and derivatives 
thereof. As mentioned in section [1 . 1 . 1 1 ChP T is an expansion in momenta (and quark masses). 
Momenta are generated by derivatives in the Lagrangian, thus we have a chain of terms with 
an increasing number of derivatives. In order to get Lorentz invariant terms there must be 
an even number of derivatives, thus we have 

£ efr = £o + £ 2 + £ 4 + - J (1-33) 

where the subscript denotes the number of derivatives. The only possible building block in 
Co is UU^ = 1 which contributes only a constant and can therefore be dropped. The lowest 
order is therefore £2 and reads 

£2 = ^(d^tf), (1.34) 
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1.4. Construction of the Effective Lagrangian 



where (...) denotes the trace which is necessary in order to get an invariant. The prefactor is 
chosen such that we recover the standard form of the kinetic term \d^(f> a d^(j) a when expanding 
U in powers of 0. The invariance of £2 under global SU(N)l x SU(N)r transformations is 
easily shown using (|l,29p and the cyclicity of traces. 

A term proportional to ((dud^U)!/') contains also two derivatives and satisfies all con- 
ditions of symmetry but is equivalent to that already written as can be shown by partial 
integration. 

Using a parametrization of the SU(N) matrices R and L we find via Noether's theorem 
the conserved vector and axial-vector currents. Expanding the axial-vector currents A^ ,a in 
the Goldstone fields we find 

A^( x ) = -F d^ a (x) + O(0 3 ), (1.35) 

which leads to a non-vanishing matrix element when evaluated between the vacuum and a 
one-Goldstone boson state: 

{Q\A^ a {x)\cj)\p)) =ip tM F e- ipx 5 ab . (1.36) 

Comparing this with (jl.22|) we conclude that Fq (the constant introduced in (|1.30p ) is the 
pion decay constant in the chiral limit. It can be measured in the pion decay tt + — > fi^Un 
(see [5]) and is found to be 92.4 MeV (actually what is measured there is F n = FQ{l + 0{m q )) 
which slightly differs from Fq since we are not living in a chiral limit world). 

Up to now we have constructed the effective Lagrangian under the assumption of perfect 
chiral symmetry. But as mentioned in section 11.3.31 chiral symmetry is explicitly broken by 
finite quark masses which mix left- and right-handed quark fields. The mass term in (jl.ip 
can be written as 

Cm = -qnMq L - q L M ] qR, (1-37) 

with M = diag(m u , m^, m s ). Although A4 is just a constant matrix, Cm would be chirally 
invariant if A4 transformed as 

M -¥ RML ] (1.38) 

(we formally replace the mass matrix with a so called spurion). One now constructs terms of 
the effective Lagrangian being invariant under (jl.29p and (jl.38p . The simplest non-constant 
term reads 

C f = El^(Mtf + UM}), (1.39) 

where Bq is an new parameter. Setting Ai now back to the constant mass matrix, the chiral 
symmetry of the effective Lagrangian is explicitly broken in the same manner as in QCD. 
The parameter Bq appearing in (|1.39|) is related to the scalar quark condensate (qq) bjQ 

F^B Q = -(uu) 3 for £17(3) (1.40) 
F 2 B = -(uu) 2 for 517(2). (1.41) 

Using (jl.39p together with (jl.34p . the masses (at lowest order) of the Goldstone bosons 
are found to be 

Ml = 2B Q m (1.42) 
M 2 K = B {m + m s ) (1.43) 

Ml = ^B (m + 2m s ), (1.44) 



2 In the case of SU(2) one writes F and B instead of Fo and Bq. 
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where m u = = m (isospin limit for simplicity) . These relations and the on-shell condition 
p 2 = M 2 are the reason why terms in the effective Lagrangian containing one quark mass 
contribute to the same order as terms with two derivatives and why therefore the effective 
mass term in (|1.39|) got the subscript 2. 

1.5 Higher Orders and Loops 

Using the same principles as in the previous section, effective Lagrangian terms of higher 
orders are constructed. As indicated in (jl.33p the next order is £4 which contains four 
derivatives, two derivatives and one power of the quark mass matrix or two powers of the 
quark mass matrix. It was written down by Gasser and Leutwyler, see [5], and is found in 
section 11.71 

Calculating now a physical process to a given order of momentum, we have to know 
which diagrams from which orders of the Lagrangian must be taken into account. This is 
what Weinberg's power counting scheme tells us. 

1.5.1 Weinberg's Power Counting Scheme 

Given an arbitrary loop diagram with L loops, I internal lines and Vd vertices of order d (i.e. 
originating from d), the amplitude will be of the form f(d 4 p) L /(p 2 ) 1 ■ Y[d(P d ) ■ The so 
called chiral order D (or chiral dimension) of the diagram counts the dimension of momenta 
and is therefore given by D = 4L — 21 + dVd- Using the identity L = I — J2d Vd + 1 one 
can eliminate / and finds 

D = J2 V d(d-2) + 2L + 2 (1.45) 

d 

This is equivalent to the following procedure: Consider a given diagram and rescale all 
external momenta as pi — > tpi and the quark masses as m q — > t 2 m q (which corresponds to 
a rescaling of the Goldstone boson masses as M — > tM). The chiral dimension D of the 
diagram with amplitude A4(pi,m q ) is then defined by M.(tpi,t 2 m q ) = t D A4(pi,m q ), see [S]. 

For small momenta only diagrams of low chiral order D will dominate. The most simple 
diagram just contains a vertex from £2 and is therefore of chiral order D = 2, i.e. 0(p 2 ), 
which is the leading order. Going on to D = 4 we will have to consider one-loop graphs 
composed only of £2 vertices and tree graphs with one £4 vertex. Thus, for fixed D there is 
only a finite number of diagrams which have to be taken into account. 

1.5.2 Renormalization Scheme 

When doing calculations at chiral order D = 4 we get loop integrals with £2 vertices which 
diverge. Since these infinities are of 0(p 4 ), they cannot be absorbed by a renormalization 
of the low energy constants Fq and Bq of £2. However, it is possible to absorb them in the 
coefficients of £4 by renormalizing just these. 

Hence, in order to get finite results at 0(p 4 ) we must not only consider one-loop graphs 
with £2 vertices but also include tree graphs obtained from £4, i.e. all possible diagrams of 
that order. There is always only a finite number of necessary 'counterterms' to cancel the 
divergences. This is meant by the statement that ChPT is renormalizable order by order. 

1.6 Local Invariance and External Fields 

So far we have considered effective Lagrangians being invariant under global chiral trans- 
formations (jl.29p only. However, in [8] it was shown that global symmetry does not suffice 
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1.6. Local Invariance and External Fields 



to determine the full low energy structure and that one needs to consider off-shell Green 
functions (i.e. involving momenta with no correspondence to the Goldstone boson masses), 
where a study of Ward-Takahashi identities (relations between Green functions originating 
from symmetries) becomes necessary. All possible Green functions may be obtained via a 
generating functional depending on external fields and the Ward identities are equivalent to 
the invariance of the generating functional under local transformations. 

Following Gasser and Leutwyler, the QCD Lagrangian is extended by coupling the quarks 
to external Hermitian matrix fields (3x3 matrices in flavor space) v^(x), a^(x), s(x) and 
p[x): 

£ = £qcd + £ cx t = £q C d + <n M (^ + 75<v)<? - q{s - h5P)q- (1-46) 

The ordinary three flavor QCD Lagrangian is recovered by setting = = p = and 
s = diag(m u , m^, m s ). The generating functional Z[v,a,s,p] mentioned above is defined as 



exp(iZ[v, a, s,p\) = (0|Texp 



i / d x£ ext (x) 



|0) (1.47) 



and is related to the vacuum-to-vacuum transition amplitude in the presence of external 
fields. Green functions are calculated by performing functional derivatives of the expression 
(|l,47p with respect to the external fields. Since the generating functional may be represented 
by the path integral 



exp(iZ[v,a, s,p\) = J[DU]exp 



d xC cS 



(1.48) 



and since we request Z[v,a,s,p] to be locally invariant, we will have to promote the global 
SU(N)l x SU(N)r symmetry of £ e ff to a local one and thus include the same external fields 
as in QCD (see [3]). 

There is another reason why it is very convenient to introduce external fields: In this way 
the coupling of an external photon field to the quarks (i.e. electromagnetic interaction) 
as well as the coupling of the massive charged weak boson Wn to the quarks (i.e. weak 
interaction) may be easily incorporated. To do so, the Lagrangian (|1.46p is split up in left- 
and right-handed parts, yielding 

£ = ^QCD + <?L7%<?L + qRl^r^qR - q R (s + ip)q L - q L {s - ip)q R , (1.49) 

where = + and L = — a^. This Lagrangian remains invariant under local 
transformation of the qi,R if the external fields are subject to transform in a certain manner, 
see e.g. [5]. 

Then, e.g. for the electromagnetic interaction we set 

V = In = -eQA^, (1.50) 

where Q = diag(2/3, — 1/3, — 1/3) is the quark charge matrix (three flavors). The weak 
interaction is included in a similar easy way. 

1.6.1 Locally Invariant Effective Lagrangian 

As in the case of gauge theories, the external fields are contained in a covariant derivative 

D IM U = d fM U-ir lx U + iUl lx , (1.51) 
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which, under a local transformation U — > R{x)UL^(x), transforms in the same way as U. 
Beside this new building block there are further objects to be used to construct a locally 
invariant effective Lagrangian: The fields strength tensors 

fjl = dp?v - <V M - i[r^ r v ] (1.52) 
f^u = d t Mlu-d u l lx -i[l IJ ,,l u ], (1.53) 

as well as the linear combination 

X = 2B (s + ip). (1.54) 

As before is of 0(p). The external fields and L are of 0{p) to match d^U; therefore 
the fields strength tensors are of 0(p 2 ) and will not contribute to the lowest order effective 
Lagrangian £2 since the Lorentz indices have to be contracted. Finally, x is °f C(p 2 ) because 

of (Q2i-(frm 

The most general, locally invariant, effective Lagrangian at lowest order now reads 

£ 2 = ^f(D,U(D^) + ^(tf/t + U X ] ). (1-55) 

It still has two free parameters; Bq is hidden in x- (|l-55p reduces to (|1.34p and (|1.39p if all 
the external fields except for s are switched off and s is set to s = M. = diag(m u , m^, m s ). 

1.7 The Chiral Lagrangian at Order 0(p A ) 

As already mentioned in section 11.51 £4 can be constructed applying the same ideas that 
where used to find £2. It does not only contain two low-energy constants but 12. Using the 
same covariant derivative and field strength tensors as above, £4 reads 

£4 = L^D^U (ZW) f > 2 + L^DJJ {D U U) ] ){D^U (D»U) ] ) 

+ L 3 (D^U {p»Trf D U U {D»U) ] ) + L 4 {D„U (D»U) ] )( X U ] + U X ] ) 
+ L 5 (D„U (D»U? {xU ] + U X ] )) + Le( X tf + U X ] ) 2 

+ L 7 ( X U^ - U X ] ) 2 + Lz{Ux ] U X ] + XU ] XU ] ) 

- iL»{f*D»U (D"U)i + {D»uf D V U) + L 10 (Uf%U^f^) 

+ Biiffrfg + f'JH + H 2 {XX ] ) (1-56) 

In the case of SU(2) the Cayley-Hamilton theorem (see e.g. [13J) can be applied to reduce 
the number of terms. 

1.8 Chiral Perturbation Theory with Baryons 

So far we have considered an effective description for Goldstone bosons only. However, it 
is possible to extend ChPT to include also nucleons (SU(2)) or the baryon octet (SU(3)). 
There is one important difficulty: The nucleon mass is a heavy mass scale that does not 
vanish in the chiral limit. The 3-momenta of the baryons need to be small in order to keep 
ChPT valid. 

The task is to find an effective Lagrangian C„n that describes the interaction between 
baryons and Goldstone bosons. We only consider matrix elements with a single baryon in 
the initial and final states. 
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1.8. Chiral Perturbation Theory with Baryons 



1.8.1 Representation of the Goldstone Bosons and the Baryons 

As in the case with Goldstone bosons we first need a suitable representation of the particles of 
interest. It turns out to be advantageous not to represent the Goldstone bosons by the familiar 
matrix U but rather by its square root u, thus u 2 = U. Under G = SU(N)l x SU(N)r this 
field transforms as 



u 



VLUrt = LuK^L, R, U) = K(L, U, R)uR\ (1.57) 



where we introduce the compensator field K(L, U, R) which is a SU(N) matrix that nontriv- 
ially depends on L, R and U, see [T7j. Only for SU(N)y transformations where L = R we 
have the simple relation K = L = R. 

We will only consider two quark flavors, thus SU(2), and therefore only the proton and 
the neutron. These can be represented by an isospinor 

-(:)■ < L58) 

The transformation behavior of the pair (U, 1 I / ) under SU(2)l x SU{2)r x U(l)y reads 

'U\ ( RUtf 



V) ' \exp{ie)K[L,R,U}^J ' (L59) 
where parametrizes U(l)y transformations. For more details see [3l I17j. 

1.8.2 Lowest-Order Effective Baryonic Lagrangian 

As in the mesonic sector, in order to incorporate Ward identities we have to construct the most 
general effective Lagrangian coupled to external fields with local SU(2)l x SU(2)r x U(1)v 
symmetry (see [5]). Local transformations imply a covariant derivative with the usual 
property to transform in the same way as \P; D fJi ^f(x) —> exp(iO(x))K[L(x), R(x), U(x)]^(x). 
The covariant derivative is given by 

= (d^ + r^, (i.60) 

where the vector 

u'idfj, - ir^u + ii(<9 M - il^)u ] (1.61) 



r M - 2 



is the so called connection which contains the same external fields as introduced in section 
1.61 Another building block is the chiral vielbein (which is an axial vector) 

u'(d l j,-ir l j)u-u(d ti -il tl )ii< . (1.62) 

Knowing the transformation behavior of under G it follows that the most general effective 

irN Lagrangian with a single nucleon in the initial and final states needs to be of the type 

'I'O^, where O is an operator transforming as O — > KOK^ under G. 

Due to different Lorentz structure of meson and baryon fields, the chiral expansion of 

C n N contains terms of all orders of p and not only of even ones, hence we have C w n = 
Ai) A2) 

itN t ^ttN ^ ■■■■ 

The most general such Lagrangian with the smallest number of derivatives being addi- 
tionally a Hermitian Lorentz scalar and even under C, P and T reads 

4jv =*(i0-m N + y7 M 75^) (1-63) 
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The nucleon mass mjv and the axial vector coupling constant cja appear as two free param- 
eters. They are both meant to be in the chiral limit (denoted by o); m = m[l + 0(m q )], 
9A = <m[1 + @( m q)] with m = 939 MeV and qa — 1-26 (known from neutron beta decay). 
The Lagrangian (jl .63|) reduces to that of a free nucleon of mass rh^ in the case of no external 
fields and no pion fields. 

The power counting rules for the new quantities are 

*,* = 0(1), D^ = G(1), {0 - mjv)* = G(p) (1.64) 

(to be explained in [5]). Especially, the fact that the covariant derivative is not counted as 
0(p) (i.e. as a small quantity) anymore is due to the fact that the nucleon mass does not 
vanish in the chiral limit and therefore the zeroth component of the partial derivative acting on 
the nucleon field does not produce a small quantity. This fact leads to an inconvenience: We 
lose the correspondence between the loop and the chiral expansion as we had it in the mesonic 
sector. The contribution from loops is not automatically suppressed and an amplitude with 
given chiral dimension D may get contributions from diagrams with arbitrary many loops. 

1.8.3 The next Order of the Baryonic Lagrangian 

(2) . 

Without any further comments on how to derive it, C ^ is quoted here (see [18): 

4at = ci(uW + u X ] u)^ - -\( u u „) (VD»D v t> + h.c.) 

4m^ 

+ ^(t^u")*¥ - ^*-f-f Uu] *. (1.65) 

Since at 0(p 2 ) only tree graphs from C *^ must be considered, there are no divergent loop 
integrals and therefore the low energy constants q cannot contain a divergent part. They can 
be determined by comparison to some irN threshold parameters which have been measured. 
The following two sets of the coupling constants are taken from |18| : 

ci = — Q.Qmjl , C2 = l^m^r 1 , C3 = — ^Am~^ 1 , C4 = 2.0mjy L 
ci = — 0.9m]^ , C2 = 2.5m^ r 1 , C3 = — 4.2m^ 1 , C4 = 2.3mjy L 



Chapter 2 

Chiral Spirals - A First Approach 



In this chapter we construct a single-particle Hamiltonian from the fermionic part of the 
effective Lagrangian. This will be diagonalized in order to find the energy of one nucleon in 
the background of pions, leading to two eigenstates with different energies. We introduce the 
spiral configuration for pions and fix the gauge such that both the connection and the vielbein 
become constant fields. Next, in the momentum space we fill the two fermionic eigenstates 
up to their corresponding Fermi energy. To do so, the sophisticated relativistic expression 
for the energy of one eigenstate is expanded in inverse powers of the nucleon mass up to 
C(m^). We express the Fermi momenta and hence the fermionic energy density in terms of 
the particle number densities of the two eigenstates. The total energy density may then be 
written in terms of these and some spiral parameters. First in the chiral limit and then away 
from it, we find critical values for the total particle number density that need to be exceeded 
in order to allow the many-particle system to lower its energy density by changing the pionic 
configuration from a homogeneous configuration to a spiral one. 

2.1 Pion-Nucleon Effective Theory 

The leading terms of the effective Lagrangian describing a pion-nucleon system with two 
quark flavors and including explicit chiral symmetry breaking due to quark masses is given 
by dEMD, dH52D and (fL63]k 

C = ^{d li Ud»rt) + ^(Mrf + UM*) +^>[iIp-rh N + ^V75?v) *• (2- 1 ) 

We will restrict ourselves to pion field configurations U (cc) which lead to a constant connection 
and vielbein, such that the Dirac equation may be solved more easily; 

u^ix) = c M , T /t (x) = d M . (2.2) 

Furthermore we consider static pion fields which implies To = no = 0. 

2.1.1 The Single-Particle Hamiltonian 

Under these conditions the fermionic part of (|2.ip is given by (dropping the symbol o on mjv 
and g A ) 

Cf = i&'fdnV + i*7*di* - mjyM + M^y^^r. (2-3) 
The equations of motion for \P are given by the Euler-Lagrange equation 

„ dCf dCt , 
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2.2. A Spiral Configuration 



which yields 

i-fdpV + i7*di* - m N y + ytSsQ* = 0. (2.5) 

The Dirac equation, on the other hand, reads 

id ^ = HV. (2.6) 

Thus, separating (|2.5|) such that on one side of the equation we have ido^f only, we can read 
off the Hamiltonian H on the other side: 

H = m Nl ° + tV** - il°fdi - ytSSsQ. (2.7) 

We have used —idi^> = pi^> (momentum of the nucleon). For H to be Hermitian we find that 
the Ci need to be Hermitian while the eZj are antihermitian, 

Ci = Ci, a T a , idi = d^ a T a with a )a , d^ a £ E. (2.8) 

Since ^ is a doublet containing two Dirac spinors, H is an 8 x 8 matrix, i.e. by the term 
"t-a^T we actually mean m7vdiag(7°, 7 ) etc. 



2.2 A Spiral Configuration 
2.2.1 The Spiral 

As explained in section [1.4.11 the pions live in the coset space SU(2)l x SU(2)r/ SU{2)y = 
SU(2) which is isomorphic (as a manifold) to S 3 . Each SU(2) matrix U can be written in 
the form 

U = £ol + ig- r, £p € R, (2.9) 

where we require the condition 

£0+^ = 1 (2-10) 

in order to have det U = 1. The condition (|2.10p explicitly shows the relation between SU (2) 
and S 3 : U can be seen as a four-dimensional unit vector. Thus, the pions (described by U) 
can be thought to live on the three-dimensional surface of the four-dimensional unit sphere. 
We now choose the following particular form for the matrix U : 

U(x) = cosaol + isinao (cos(^(x)r 1 + sin </?(x)t 2 ) , <f(x) = a ■ x + ipo, (2-H) 

which satisfies the condition (|2.10p , a is a constant three-dimensional vector. In appendix 
IA1 it is shown that the condition of having a constant connection as well as a constant mass 
term in the Lagrangian naturally leads us to the configuration (|2,lip . 

Let us interpret this configuration: Since we have £3 = we restrict ourselves to the two- 
dimensional hypersurface Co + ^1 + ^2 = 1> i- e - to S 2 . Hence, the matrix describing the pion 
fields in each point of the three-dimensional space can be respresented by a three-dimensional 
vector lying on the unit sphere at that point. The ground state (no pions) corresponds to 
U = 1 (see section II ,4.2[) which is equivalent to ao = 0. So, ao may be interpreted as the 
polar angle of the unit sphere, if we choose the ground state to point in the third direction 
of a cartesian coordinate system. Then, <p(x) can be thought of as the azimuthal angle and 
we can indeed describe the vector representing U by the ordinary spherical coordinates. 

Now let us walk around in space, namely in a plane perpendicular to the vector a. This 
plane may be parametrized as x = xo + Xv + fiw where xq, v and w are constant three vectors 
(v and w linearly independent) and A, /i £ K. Since a is perpendicular to both v and w we 
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get x ■ a = xo ■ a = const. Thus, in the plane perpendicular to a the azimuth angle if (see 
()2.1ip ) does not depend on x and the unit vector respresenting U will not change throughout 
space. However, if we walk in any other direction the azimuth angle f{x) becomes position- 
dependent and the vector starts to rotate about the "north pole" of the unit sphere. For 
a given speed of our walk the scalar product a ■ x and so the azimuth angle increases most 
quickly in the direction of a which we therefore may identify as the direction of the spiral. 

Having such a configuration, we can choose an arbitrary vector to represent the ground 
state and hence to define our third direction. In all neighboring space points (except those 
lying on the plane perpendicular to a) the vector will have rotated. 

For our calculations we need the field u rather than U . We make the ansatz 

u(x) = cos at + i sin a e(x) ■ f , (2-12) 

where e(x) is a unit vector. Using (e • r) 2 = ejej-rV- 7 = eiej(ie^ k T k + 5^1) = ejejl = 1 
(where we used the symmetry of and the antisymmetry of e ijfe in the indices i and j and 
the fact that e is a unit vector) we find 

u 2 (x) = (cos 2 a — sin 2 a) 1 + 2i sin a cos a e(x) ■ r = cos 2at + i sin 2a e(x) ■ r 

and the requirement u 2 = U yields 

a = —. e{x) = (cos tp(x), sin ip(x), 0). (2.13) 

2.2.2 Connection and Vielbein 

From this we may calculate the connection and the vielbein (we set the external fields and 
to zero) . For convenience we write = 5u cos <p + b2i sin <p and, using the particular form 
of if of (|2.1ip . diej = ai(S2j cos ip — 5\j sin ip). We then have 

= - u^diU + udiv) 

= — [(cos al — i sin a eiT l )i sin adiCjT^ — (cos at + i sin a eiT % )i sin adiejT^ 

= — ■ 2 sin 2 aekT k diejT^ = en sin 2 a(5\k cos if + 82k sin f) (faj cos f — 8\j sin if)r k T^ 

= en sin 2 a [^ifc^j cos 2 f — 5\k8ij cos ip sin f + b2k&2j sin if cos ip 

-5 2 kSi jS m 2 f](ie k ^T l + 5 k H). 

The second and third term in the squared bracket vanish because of (anti-)symmetry of the 
indices k and j when multiplied with e k ^ 1 and because they cancel when multiplied with 5 k i . 
The first and the fourth term disappear when multiplied with 8 i because the indices can not 
coincide in both Kronecker deltas at the same time. Thus, what is left is 

Ti = iai sin 2 q(cos 2 fe l2l r l - sin 2 fe 21l r l ) = im sin 2 err 3 . (2.14) 

The vielbein: 

v)diU — udiiJ 

[(cosal — i sin a ei r l )i sin adiejT 3 + (cosal + i sin a e{T l )i sin adiejT J ~^ 
-2 cos a sin adiejT^ = — 2aj cos a sina(— sin fT 1 + cos ipr 2 ). (2-15) 



Ui = i 
= i 
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Obviously U{ is not a constant field in contrast to what we demanded in (|2.2|) . It can however 
be made constant by an appropriate SU(2)y gauge transformation. As can be shown, the 
transformation behavior of it M under local SU(2)l x SU(2)r is given by 

Up^Kuprt, (2.16) 

where K is the compensator field introduced in section 11.8.11 Starting from the equation 

d^K = KT^ - T'^K (2.17) 

(derived in [5]) and using the relation d^KK^ = —Kd^K^, we find the transformation be- 
havior of the connection to be 

r M Kpp + dJK*. (2.18) 

In the case of the isospin transformation L = R = V we have K = V (see |5j). Writing out 
Uj in matrix form, 



Ui{x) = —2<ii cos a sin a 
= —2a,' cos ol sin o. 



— sin</?(x) — i cos ip(x] 

s\mp{x) + icosip{x) 

— i exp(— iip(x)) 

iexp(i(p(x)) 



and parametrizing the SU (2) matrix V in the following way, 

v(x) - (eMix(x)) \ 

we find, using (j2.16fl . for the transformed vielbein 

if \ n ( —iexp(—i<p(x) + 2ix(x))\ ,„ on > 

U ^ = - 2aiCOSaSma {ieMi<p(*)-2ixW J (2 - 20) 

and, using (j2. 18fl . for the transformed connection 

Ti(x) = iai sin 2 ar 3 - idix(x)r 3 . (2.21) 

We have dropped the arrow on x and will go on like this for convenience (remember that we 
do not have time evolution). (|2.20p leads to the choice x( x ) = ^f( x )\ then the vielbein gets 
constant and we obtain 

= mj(sin 2 a — ij) T ^ = d% = const (2.22) 
u\ = —2di cos a sin ar 2 = a = const. (2.23) 

For (|2.22p we have used the form of (p(x) introduced in (|2.1ip . 
2.2.3 Pion Contribution to the Energy Density 

Making use of U = u 2 , v)u = 1, the identitie^j] ud^v) = —d^uv) and u^d^u = —d^vf 
well as the cyclicity of the trace, it is easy to show that 



•ull U as 



{dffldPrf) = {u^). (2.24) 



1 both following from uu* = u'u = 1 
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With ([2,23|) we thus find (dropping the prime) 

{d^Ud^) = {uy) = 4 cos 2 a sin 2 cm; ^(1) = 8 cos 2 a sin 2 aa^ . (2.25) 
We finally need 

{Mrf + UM j ) = {M{U j + U)) = 2 cos «o {-M. ) = 2 cos ao (m u + rrid) , (2.26) 

where we have used JVv = Ai. Thus, remembering (|2.13p . the mesonic part of the Lagrangian 
reads 

C p = 2F 2 cos 2 a sin 2 actio 1 + F 2 B cos ao{m u + ma) 
F 2 

= — sin 2 aoa-ia 1 + F 2 B cos ao(m u + m^). (2.27) 

Since we have 8qU = we obtain for the energy density due to the pions 

F 2 

£p = Hp = —C p = — — sin 2 a^cua 1 — F 2 B cos ao(m u + m^) 
F 2 

= -^o 2 s™ 2 a o — F 2 B cos ao(m u + m^), (2.28) 
where in the last step we used OjO* = —a^ai = — |o| 2 =: —a 2 . 

2.3 Fermionic Contribution to the Energy 

The fermionic energy is found by calculating the eigenvalues of the single-particle Hamiltonian 
(|2.7p . In principle this corresponds to the calculation done in section 12.2.31 for pions but since 
these are scalar degrees of freedom, this calculation was much simpler. 

Comparing (|2.8|) with (|2.22p and (|2.23j) . we find that only a t 2 / and d^z / and we 

set 

Pi '■= Q,2 = — 2ajCosasina = — Ojsinao (2.29) 

Si ■= d i;3 = -Oj(sin 2 a - -) = y cos a Q . (2.30) 

The Hamiltonian can then be written as 

H = m Nl °t + 7 Vpil " l°7%r 3 - ^tVtsA^ 2 , (2.31) 

where 1 means the 2x2 unit matrix. The eigenvalues of the Hamiltonian are found by 
Mathematica and read 



Et = \Jp 2 + di + m 2 N + l -g\^ ± ^{pA? + g 2 A {m 2 N [3 2 + & 2 P 2 + faft) 2 - (W 2 } (2.32) 

Eg = -^P 2 + 5* + m 2 N + \g\^ ± y^^) 2 + g\ {m 2 N ^ + + foft) 2 ~ ( W 2 }, 

(2.33) 

where p 2 = p\ + p\ + p 2 etc. Each of these four eigenvalues is twofold degenerate, such that 
we have eight eigenvalues altogether. Although in the original basis \& has the form (jl.58|) . 
in the new basis where H is diagonal, the eigenvectors contain linear combinations of the 
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different p- and n spinor components. It is therefore not straightforward to decide which 
eigenvalue belongs to which state. However, negative energy states (which are the antistates 
to the positive energy states) should not be considered in the effective theory, since we work 
at an energy scale of and consider only small deviations from it. Considering the negative 
energies would mean a deviation of 2m jv which is way to large to do ChPT. Hence we neglect 
the E±. 

Due to (|2.29p and (|2.3U|) the two vectors /3 and 5 are parallel and therefore the two terms 
5 2 (3 2 and (5if3i) 2 = (5 • j3) 2 cancel, such that the fermionic energy reads 



E± = ^p* + 6 2 +m 2 N + ^2 ± ^4^. §) 2 + g2 A { m 2 N 02 + tf. £)2 } . (2 . 34 ) 

Making an expansion of (|2.34p in powers of 1/mjy yields 

V 2 5 2 1 / 1 \ 

E± = m N + -f— + ± -g A /3 + O — , (2.35) 

2rriN 2totv 2 V m w / 



where /3 = \J 0\ + /3| + /?f > 0. S-t is the energy of one linear combined state. We now 
want to calculate the energy due to a given amount of states. To do this, we fill, in a given 
volume L 3 , the fermions up to some Fermi energy. Let us suppose the temperature T = or 
at least a sufficiently low temperature to assume that we have occupied states only up to the 
Fermi energy. Since the dispersion relation (|2.35|) depends only on p 2 , in momentum space 
the Fermi surface is a sphere and to the Fermi energy corresponds a Fermi momentum pp. 
The number of possible modes in a region d?p is 



2ir J 

which (neglecting of course the interaction between the fermions) remains the same even if we 
consider two different particles (since all the creation and annihilation operators of different 
particles anticommute with each other, the statistic of one particle is completely independent 
of the other). Hence, we can talk about two separate Fermi momenta pf+ and pf-- Given 
these, for the particle numbers N± we then find 

Fermi 
sphere 

The factor 2 appears because the states E± are each twofold degenerate. We calculate N± 
since we want to express the fermionic energy in terms of the particle number or particle 
number density respectively rather than in terms of Fermi momenta. The particle number 
density is 

N± Pf± ,„ ^ 

= TF = 37T 2 "' (2 - 37) 



Now, using (|2.35p the total energy of all fermions of a given sort up to the Fermi momentum 
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is 

f L \ 



Fermi 
sphere 



m v I -- ± -fl, M 2 fi-V f d * p+2 (lL) 3 [d 3 p 1,2 



2m j\r 2 / \2ir J J \2tt J J 2niN 

Fermi Fermi 
sphere sphere 



m - + 2^ ± \ UP ) N± + ^~ N it) ^ I dPp4 

o 

r2 -1 \ 7~ 3 p 5 
m N + - ± -g A N± + — 2 (2.38) 

Using now (|2.37p to replace pp± in terms of n± we find for the energy densities 

+ EL ( 52 1 n\ (3vr 2 n±) 5 / 3 , . 

The total fermionic contribution to the energy density e t = ef ot + e^ t finally reads 

/ 5 2 1 \ (3^ 2 n+) 5 / 3 / 5 2 1 \ (3vr 2 n_) 5 / 3 

e/ = mjv + + 5A /3 n + —— + [m N + - -g A p n_ + — — 

\ 2rri]y 2 J lUir z mN \ 2ra^ 2 ) Vdix^m^ 

, 5 2 \ g A p (3^ 2 ) 5 / 3 / 5/3 , 5/3\ (i) An , 

mN + ^) n + — {n+ - n - ) + m^\ n ^ +nj ')> (2 ' 40) 

where in the last step we set the total particle number density n = n + + n_ . 

2.4 The Total Energy Density 

The total energy density is given by this fermionic part and the pionic part of (|2.28|) . As- 
suming for simplicity the isospin limit m u = we can use (jl.42p . Together with (|2,29p and 
(|2,30|) and remembering that n + + n_ = n fixed and so n_ = n — n+ we then get for the 
total energy density the following expression: 

,2 „„„2 



/ a cos ao\ g A 
etot = e/ + (-p = \ rn N H — In- — asma (n - 2n + ) 



L 2 

W7r 2 m N Y" + ' v " ) ' 2 

Let us find the value of a that minimizes e to t : 



+ (»f + (»-» + ) 5/3 ) + sin 2 n„ - f' 2 3/3<osn„. (2. !J ) 



<9e to t an cos a . , x , ^2 • 2 ! n 

— sin ao(n — 2n + ) + ^ asm ao = U 



= 9 r~TZ ^9 — 2 ' ( 2 ' 42 ) 



9a 4m jy 2 

2ttin g A sinao(n — 2n H 
n cos 2 ao + 4m-ArF 2 sin 2 ao 

The second derivative is always positive for < ao < vr/2, showing that we have found a 
minimum. After plugging a m \ a into e to t the first derivative with respect to ao reads 

<9e to t ^2^2 • 2o 2 1 mA r n(n-2n + ) 2 sin2ao , 

— — = F M_ sin ao — — — ^ -77. (2.43 

<9a (4F 2 m7v + n + (n-4F 2 mAr)cos2a ) 2 
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For a minimum this expression needs to vanish and we get 

F 2 M 2 sina (AF 2 m N +n+(n- AF 2 m N ) cos2a ) 2 = 2g\m N n(n - 2n + ) 2 sin2a . (2.44) 
This equation determines how the angle ao will adjust itself for given densities n and n + . 

2.5 Chiral Limit 

Let us now consider the chiral limit, i.e. set M n = 0. Then the left hand side of (|2,44p 
vanishes and for n + ^ re/2 we conclud^l sin2«o = 0. In the interval < ao < vr/2 this has 
two solutions; ao = and ao = 7r/2. Let us examine the second derivative of e tot with respect 
to ao at ao = vr/2: 

d 2 € to t 



8oq 



= g 2 A n(n - 2n+) 2 
a =^ 16F 4 m/v V ' ' 



and 



This is always positive (except for n + = re/2 where it vanishes) and thus the minimum of 
etot is at ao = vr/2. So, in the chiral limit the vector representing the pions will move in the 
equatorial plane. There (|2.29[) and (|2.30p become 5i = and = — aj, hence a 2 = j3 2 . We 
therefore get 

e P = ^a 2 = —P 2 (2.46) 

e f = m N n - 9 -^-(n - 2n + ) + (ref + (re - n+) 5 / 3 ) . (2.47) 

Let us find the value for /? such that etot = e y> + e / gets minimized: 

^ = F 2 ?-f (re-2re + ) = => P = §- 2 (n-2n + ) (2.48) 

Since f3 > we must have re+ < n/2, which makes sense as the lower energy states EL are 
energetically favored and n + > n/2 would only increase the energy density. Plugging the 
expression for (3 into e to t yields 

e to t =m N n- 9 ^(n- 2n + ) 2 + g£ (ref + (re - re + )^) . (2 . 49 ) 
At a fixed total density re we now vary re+ in order to find extrema of the total energy density: 
detot (3tt 2 ) 5 / 3 / 2/3 , 9a^ x ! n fn ^\ 

From a plot of det t/dn + (see fig. I2.1a|) one finds that for small re this equation has the only 
solution re + = n_ = re/2. The second derivative at re+ = re/2 reads 



d 2 e to t 



drij 



97T 2 m]\f V re / F 2 



which is positive for 



16F 6 vr 

n<n c = — ^ — s-. (2.52) 
3g A m N 



2 For n+ = n/2 we have a = according to l|2.42|) . However, this implies a non-spiral configuration, see the 
remark preceding (|2,53[) . 
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The symbol o reminds us that we are in the chiral limit. 

Thus, for densities n smaller than this value n + = n_ = n/2 is the only minimum of etot- 
In this case we have, using (|2.48j) . (3 = 0. But since f3 = a = \a\ it follows that a = and so 
ip{x) = ipo (see (|2.1ip ). In this case we therefore do not have a spiral. The energy density of 
this non-spiral configuration is 



£ » = '"""+ 5^257^=™"" +1^(4'") ■ 





Figure 2.1: In the upper line detot/dn+ is plotted against n + /n. The lower line shows 
e* ot (n + ) = — ^^+=1/2) ~ ^' mu ltip ne d by an appropriate constant indicated in the plot. The 
horizontal line in the lower plots marks ej^ t and shows that for n > h c a spiral configuration 
is favored energetically. 



However, for densities n > h c an interesting thing happens: Since in this case the second 
derivative of etot at n + = n/2 becomes negative, n + = n/2 becomes a maximum. Addition- 
ally, (|2.50p now has two roots (see fig. I2.1bp . where, since the one at n+ = n/2 is a maximum, 
the other is a minimum. 

With increasing n the additional root of (|2.50p moves to the left of the interval < n + < 
n/2 and eventually leaves it. The "upper critical density" where this happens is found by 
solving detot /dn+ (0) = for n, which yields 

hi = V^T- (2-54) 
9A m N 

For densities n > n" det t/dn + looks like in fig. I2.1cl and we do not have a local minimum 
in the interval < n + < n/2 anymore but instead a global one at n + = 0. Hence, is the 
density from where on E_ is completely filled. 

The relevant fact is that for n > h c the energy density gets lower than for some 
n + 7^ n/2 and hence for a spiral configuration. 

For the calculation of the numerical value of h c we use, as mentioned in the introduction, 
the values of the low energy constants in the chiral limit, as well as the physical ones. The 
value of the pion decay constant in the chiral limit was found in [14] . we take F = 86MeV, 
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2.6. Away from the Chiral Limit 



while the physical value is taken to be F = 93MeV. The nucleon mass in the chiral 
limit was calculated in |19j : we set rh^ = 882 MeV. For the physical nucleon mass we choose 
rriN = 939 MeV (neutron mass). Finally, the axial- vector coupling constant §a in the chiral 
limit was extrapolated in [20J with a rather large uncertainty. We take §a = 1-20, while the 
physical value is more precisely known; we take qa = 1-26. 

If h c denotes the critical particle number density in the chiral limit, (|2.52p then yields 



102.6 GeV 3 or n 1 / 3 = 46.81 MeV (2.55) 
when using the values of the low energy constants in the chiral limit and 

101.5 GeV 3 or h}J 3 = 46.64 MeV (2.56) 



when calculating with the physical low energy constants. We see that the difference between 
these two results is quite small, i.e. the value of the critical particle density does not strongly 
depend on uncertainties in the low energy constants. This will also be true for all the further 
results and for this reason we will quote the result obtained using the low energy constants 
in the chiral limit only. 

Let us compare the value in (I2.55P to the nuclear matter density. The radius of a typical 
nucleus of an atom with A nucleons is R ~ A 1 / 3 tq 1 where tq ~ 1.3 fm. With 



4-7T n 

A = —nR 3 
3 



it follows for the nuclear matter density 



3 

"-nuclear = ^3 = 0.11 fnT 3 , (2.57) 

or, using 1 f m = 1/197.3 MeV -1 , 

^udear « 95 MeV. (2.58) 

Hence, at least in the chiral limit, we can expect a spiral configuration to occur already at 
densities smaller than ordinary nuclear matter. 



2.6 Away from the Chiral Limit 

As soon as we allow the quarks to bear masses, M n will no longer be zero and the vector 
representing the pions will no longer move in the equatorial plane. Rather, ceo will adjust 
according to (|2.44|) . 

Plugging ()2.42p into ()2.4ip and varying the resulting expression with respect to n + leads 
to the equation 

5(9 7 r 4 ) 1 / 3 ncos 2 a (nf - (n - n + ) 2 ' 3 ) 

+ 4m 7 vsin 2 ao(c/im i v(n-2n + ) + (97r 4 ) 1/3 F 2 ^n+ /3 -(n-n + ) 2/3 ^ =0. (2.59) 

Again, this equation is always solved by n + = n/2. According to ()2.42p this means a = 
and hence a non-spiral configuration. 

Since we have two parameters now, n + and ao, there does not seem to exist a straight- 
forward way to analytically determine a critical fermionic density from where on a spiral 
configuration (n + ^ n/2) is energetically favored. Hence, we approach this problem numeri- 
cally. 
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The energy density of reference we aim to undercut is the one of the non-spiral-configuration 
where n+ = n/2 and a = 0. From ()2.4ip we find 

3 f§ \ 1 / 3 

e+otl n ln =rriMn-\ ( -7r 4 n 5 ) — F 2 M^cosan, (2.60) 

tot la=0,n+=n/2 ^ \®m N \A J ' V ' 

which is minimized by ao = 0. The non-spiral energy density therefore reads 

1 /3 

If we take the total energy density as a function of the two parameters, €tot(n+, cto), then the 
interesting quantity to look at is 

d(n + ,a ) := e^ t - e tot (n + ,a ) = e to t(n/2,0) - e to t(n + ,a ), (2.62) 

which denotes the difference of energy densities between the non-spiral configuration and 
a particular configuration with parameters n+ and (Xq. Whenever d(n + ,ao) < 0, the non- 
spiral configuration will be energetically favored. But as soon as we can find a configuration 
such that d(n + ,ao) > 0, this configuration will have n + ^ n/2, since (|2.60p tells us that 
for n + = n/2 the parameter ao will go to zero and hence d = ^ 0. This configuration 
will therefore be a spiral configuration and will have a lower total energy density than the 
non-spiral one. The fermionic density where this is possible for exactly one configuration will 
be our critical density n c . Let us refer to the corresponding parameters as the critical ones; 
n% and a%. Fig. O shows the function d(n + ,a ) for M n = 2MeV and n 1 / 3 = 49.7MeV, 
using the low energy constants in the chiral limit. 




Figure 2.2: d(n + ,a ) for M n = 2MeV and n 1 / 3 = 49.7MeV 

Obviously we are in this case already above the critical density since d(n + ,ao) > for a 
whole region in the parameter space. 

A first observation is the following: For small values of the critical parameter n c + will 
lie somewhere in the interval < < n/2. For increasing values it will move towards lower 
values and eventually reach n\ = for a certain pion mass M%. This is visualized in fig. 
12.31 where contour plots of d(n + , ao) for different values of M T just above the critical density 
n c are shown. For comparison: Fig. I2.3al shows the same situation as fig. 12.21 From this 
plots the above statement is obvious. For pion masses larger than M£ ~ 20MeV the critical 
parameter n c + drops down to zero. For an investigation including physical pion masses we 
can therefore clearly constrain our analysis to the n + = parameter subspace. 

Although we now have just one parameter ao to consider, the problem of finding the 
critical fermionic density n c from where on et o t(0, ao) < £tot( re /2, 0) still seems to be nontrivial 
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Figure 2.3: Contour plots of d(n_|_,ao) ° ver the parameter space for different pion masses. 
The fermionic densities are chosen to be slightly above the critical ones. 



(since the condition <9e to t(0, oto)/dao = is complicated to solve for ckq). We therefore go 
on numerically in the following way: For each value of n in a given interval we numerically 
calculate the maximum value of d(0, (Xq), which we call c2 max (n). As long as d max < 0, the 
non-spiral configuration will have a lower energy density. For large enough values of n, we 
will however find d max > 0, which means, that the spiral configuration will appear. The 
critical density n c is therefore just the root of d mscK (n). 

The plot in fig. I2.4al shows d max (n), numerically found by Mathematica, using a physical 
pion mass of M n = 135 MeV. The continuous graph was obtained by using the low energy 
constants in the chiral limit, while the dashed graph shows d max (n) using the physical values 
(again showing the weak dependence of the result on variations of the low energy constants). 
For small values of n d max is small but negative, actually d max — > as n — > 0. This is 
consistent, since at small densities the system will increase its total energy density only 
slightly when choosing a spiral configuration instead of a homogeneous one. For increasing 
particle density this effect however becomes stronger, i.e. the difference between etot an d e^ t 
becomes larger. However, for densities larger than around 85 3 MeV 3 the difference starts to 
shrink again since the lowering of the fermionic energy due to the spiral configuration (that 
increases the energy in the pionic sector) starts to carry weight. For densities higher than 
around 100 3 MeV 3 the decrease of the fermionic energy becomes larger than the increase of 
the pionic energy, i.e. the total energy density is lowered for a spiral configuration and d max 
becomes positive. The root, i.e. the critical density is found to be located at 

n c = 1069.6 GeV 3 or n l J 3 = 102.27 MeV. (2.63) 

Beside vanishing and physical pion mass, let us find the critical density additionally for 
some intermediate pion mass. We take M n = 40 MeV, which is considerably lower than the 
physical mass but still higher than the critical mass of about 20 MeV from where on we can 
set n + = 0. Fig. I2.4bl shows again the function <i max (n). The critical density (let us denote 
it by n c ) is 

n c = 345.4 GeV 3 or n l J 3 = 70.16 MeV. (2.64) 
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n 1 ' 3 [MeV] n 1 / 3 [MeV] 

(a) M-n = 135 MeV (b) = 40 MeV 

Figure 2.4: The maximum values of d(0, «o) as a function of n for physical pion mass as well 
as for M n = 40 MeV. Continuous graph: Low energy constants in the chiral limit; Dashed 
graph: Physical low energy constants. 



Chapter 3 

Higher Orders 



The calculations in chapter [2] showed that the critical particle number density for the spiral 
configuration to occur lies in the region of nuclear matter density. It is not obvious that we 
can still trust ChPT at such high fermionic densities. In this chapter we therefore want to 
take higher orders into account and check their influence on the results. First, we continue the 
expansion of the relativistic expression for the fermionic one-particle energy up to 0{m~^). 
This will point out the so far uncared problem that the expansion in the way we do it is not 
valid for all values of the spiral parameters. After having realized that the simple continuation 
of the energy expansion will yield no results, we perform the expansion explicitly for small 
values of the parameter ceo- 

Next, we take higher orders of the baryonic Lagrangian into account, breaking chiral 
symmetry explicitly also in this sector. 



3.1 Higher Orders in the Energy Expansion 

First, let us examine how the above results will be affected if we continue the non-relativistic 
expansion ()2.35|) of E± up to order l/m 2 N . We then have 

P 2 S 2 1 n 

E±=m N + h ± -g A /3 

Lrrtfq Lrrtfq 2 

± < 4(? - Sf + A ^ Bf - ™ + ^ } + (4) ' (3 ' 1) 

The first term proportional to shows a severe problem: The j3 in the denominator may 
get small for small values of ato and the expansion fails as the term explodes. This will be 
analyzed further, see below. 

Unlike E± of f|2.35[) the energies now not only depend on the absolute value p of the 
fermion momentum, but also on its direction. The Euclidean scalar products p- 5 = pibi and 
p ■ f3 = Pifii may be rewritten in terms of the absolute values of p, 5 and (3 and the angles 
Z(p, 8) and Z(p,f3) respectively, which turn out to be the same, namely Z(p,a), because of 
(I2.29P and (|2.30p . Denoting this angle by 6 and choosing the coordinate system such that 
the third axis points into the direction of a, we have 

p ■ 5 = p5 cos 9 = p3<5 (3-2) 
p ■ f3 = p/3 cos 6 = ps/3. (3.3) 

Since E± is now not only a function of p 2 but also of 9 or p% explicitly, the Fermi surface is 
not a sphere anymore. In order to see its shape in the momentum space, let us rewrite ()3.ip 
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as 

E± = wf + wf p 2 + wfpj, (3.4) 

where 

w =m N + ± -g A p =F -. 2 (3.5) 

W 2 = n T 1 T ( 3 - 6 ) 

2 1 „2 o2 



,±_ ^4^+^/3 



~ 3 " ^ 4^/^ (3 - 7) 

are constants with respect to pj. Fixing now -E^ to certain values (e.g. the Fermi energies) 
Ep, we may write f)3.4[) as 

W 2 D 2 , ™f p 2 , 2 _ x f o 8) 



which is the equation of an ellipsoid in momentum space. A general ellipsoid 

a? b 2 c 2 



>A + I + 4 = 1 (3.9) 



has the volume 

47T 

Ellipsoid = -^-abc, (3.10) 
so the volume of our Fermi ellipsoid is 

(3.11, 



• \Jw^wt 



When calculating the particle numbers iV-j- (cf. ()2.36p ) and the total energies E^ t (cf. (|2.38p ) 
we have to integrate over this Fermi ellipsoid. Let us first find the particle number: 

ellipsoid 

In order to calculate the total energy of all fermions of a given sort, we have to integrate E± of 
fj3.4[) over the Fermi ellipsoid. From (j3.8|) we see that the integration variable transformation 

2 > 2 
Pi ~>Pi 

2 > 2 

P2 ~>P2 

Pl^p'i = ^^pl (3.13) 
^2 

maps the ellipsoid to a sphere S of radius 
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which is an integration region we can easily deal with by changing to spherical coordinates; 
pi = p sin 6 cos <p, p2 = p sin 6 sin cj), p% = p cos 9. r$± is the analogon of pf± in section 12.31 

EL = 2 (^) 3 / d 3 pE±(p) = 2 (^) 3 / d 3 P {wf + wfp 2 + wfp 2 ) 

Fermi Fermi 
ellipsoid ellipsoid 



Fermi Fermi 
ellipsoid ellipsoid 



N± 



= tof W± + ^'4 ± w}J ± "f ± . (3.15) 

Now we replace rs-t such that we have the total energy in terms of the particle number 
(density) , in analogy to what we did in section 12. 3[ Dividing (|3.12|) by L 3 we get the particle 
number densities n± which we may solve for Ep — wf. Plugging this expression into (|3,14p . 
we find 

1/3 



r s± = ^2 1/6 i^Tr 2 n ± ■ sjwf + wf\ . (3.16) 
After using this to eliminate r$± in (|3.15|) and dividing by L 3 we are led to 

4,=»fn ± + ^4 / V( 1 + ^|) 1/3 . (3.17) 

Let us split up the terms wf; wf = xi ± ?/«. We have (see (|3.5p - (|3.7p ) 

<5 2 , , 

xi = m N + - (3.18) 

IrriN 

12 - (3 ' 19> 

x 3 = (3.20) 



W -5«^ I'-inT 1 (3 ' 21) 



IB = -|4 (3-22) 
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The total fermionic contribution may then be written as 



w^n + + w 1 n_ + 



x\n + yi(n + - nJ) + 



35/3^4/3 



5 / 3 „„+ 



nj_ w 



1 + 



Wo 



1/3 



5/3 - 
+ n_ u; 2 



(Co 



- \ 1/3' 



1 + -^ 



W 



35/3^4/3 



(x 2 + y 2 ) 1 + 



2/3 



X2 + 2/2 

+ (x2 - 2/2) ( 1 



1/3 



n 



2/3 



5/3 



X2 ~ 2/2 



1/3 



71 



5/3 



(3.24) 



Since we are only interested in terms up to order we expand 



(X 2 ± 2/2) 1 ± 



2/3 



X2 ± 2/2 



1/3 



in powers of m^ 1 : 



(X 2 ± 2/2) 1 ± 



2/3 



X2 ± 2/2 



1/3 



1 

2m7v 6m 2 



+ 



5 2 



iV 

§2/2 



+0 



in 



A' 



x 2 ± y 2 ± 



2/3 



(3.25) 



J3 1 M 
3 3 



Plugging this into f|3.24|) yields 
ef = x\n + yi(n+ - n_) 



+ 



35/3^4/3 



5 

-2 



, 5/3 . 5/3 
X2 I n+ + n_ 



+ (2/2 + y 



5/3 5/3 

n ! — n_ 



(3.26) 



When neglecting all terms ~ m N this expression reduces to that of (|2.40p . i.e. the expansion 
is consistent. 

Using (|2T29l) and (pDtOD . let us study the term y 2 + y 3 /3: 



2/3 
2/2 + — = 2/2 + ~ — o 2 _ 
3 3g A pm z N 



2/2 

3 



5A /3 _ 25 2 -g\p 2 



3g A /3m 2 N 6m 2 N 6g A /3 



m 



N 



(cos ao cot ao — 2g A sin ao) 



(3.27) 



12g A m% 

For ao — J- this term diverges to +00. Since (n ! — n_ ) < for < n + < n/2 it follows 
that €f — > —00 for ao ~~ ► 0. Obviously, this behavior is non-physical, namely the energy 
would no longer be constrained from below. Here, we are confronted with the problem we 
already mentioned at the beginning of this section, namely that the expansion of the fermionic 
energy in the way we did it is valid only for ao not too small. A more detailed analysis of 
the expansion (see appendix [B]) yields the condition 



2m^g A sinao 



< 1 



(3.28) 



for the expansion to valid. Since the parameter a will be known only after minimizing the 
energy with respect to it, it is not a priori possible to decide down to which ao we may trust 
the expansion. We therefore now first go on using (|3.26p and find out whether some relevant 
information may be extracted. 
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3.1.1 Expansion for a ~ 0(1) 

Besides ()3.26p . for the pionic contribution we use the same expression as before, i.e. (|2.28[) 
together with f| 1 . 42 j) . The total energy density is then again given by etot = e / + e p- After 
replacing the x% and yi with the expressions of (|3. 18[) - ()3.23[) the expression reads 

£ » = ( mK + sb) " + W (n + - n - ) (' - 4£) 

+ -H^l Lf + „«/» + (nT _ „ 5 /3^ 

F 2 

+ — a 2 sin 2 a - F 2 M 2 cosao. (3.29) 

Finally, we use (|2.29p and (|2.30p to replace /3 and 5 and substitute n_ = n — n + to get the 
total energy density as a function of a, ao, n and n + . As we did in section 12.41 we now 
minimize etot with respect to a. The equation detot/da = is a quadratic equation in a 
possessing the two solutions 

a± = ± 2 — — ^ — (=f20g A m N (n cos 2 a + ^F 2 m N sin 2 a ) + 5 1 ) , (3.30) 

6Ug A {n — 2n + ) sin ao cos z ao 

where S is the square root 

{/ o o o o \ 2 

400 (g^mArncos ao + 4F gA^N sm a o) 

— 240^(72 — 2n + ) cos 2 ao sinao (97r 4 ) 1//3 ( nV 3 — (n — n + ) 5//3 ) cos ao cot ao 



x 1/2 

+ 2 5 isinao{(9vr 4 ) 1 /3^ n _ n+ )5/3_ n 5/3^ _ 5m2v(n _ 2n+) jj I _ (3 31) 

Since a_ < and a+ > for all values of n, < n + < n/2 and < ao < vr/2, we conclude 
that a_|_ is the physical solution, d 2 e^t / da 2 (a+) > confirms that a + minimizes etot- Hence, 
plugging a + into etot we eliminate the parameter a and are left with three parameters. 
Next, we identify the non-spiral energy. Using de l'Hopital we find 

lim a+ = 0, (3.32) 

so as before the spiral disappears for n + = n/2. Setting n + = n/2 and then a = in (|3,29p 
we find 

3 /9 \ 1 ^ 

etotl /o „ = miviiH I -7r 4 n 5 ) — F 2 M 2 cosao, (3.33) 

tot| n+=n/2ja=0 -r 10m ^ i 4 i o, v 7 



1 /3 

e^t = + f ^vrV) - F 2 M 2 , (3.34) 



which is minimized for ao = 0. Hence we can write 

3 /9 
Wm N \4 

which is the same expression as (|2,6ip . 

The two plots in fig. I3.1al and I3.1bl show the energy densities both of lower and of higher 
order as functions of ao, as well as the non-spiral energy density (horizontal line). The particle 
density was chosen to be n = 102.27 MeV, i.e. the critical density we found in section [2T6l 
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(a) n+ = (b) n+ = O.ln 



Figure 3.1: Lower and higher order energy densities at n 1 / 3 = 102.27 MeV, at two different 
values of n+, using the physical pion mass 

The plot in fig. I3.1al shows the situation for n + = 0. The minimum of the lower order 
energy density (upper curve) touches ef^, confirming that the particle density is set to the 
the critical one of the lower order calculation. The higher order energy density (lower curve) 
however shows the problem we already predicted earlier, when having analyzed expression 
(|3.27|) : The function does not show a minimum but it rather diverges towards — oo for small 
values of ao- This unphysical behavior is due to the fact that the expansion of the fermionic 
energy in the way we did it in (|3.ip breaks down for too small angles. This is confirmed also 
quantitatively: Using (|3.30p we find that in the case of n 1 / 3 = 102.27 MeV for a < 0.13 the 
left hand side of (|3.28|) becomes larger than 1 and hence the expansion (|3.1|) starts to fail. 
This situation makes it impossible to identify a minimum of the higher order energy density 
and to decide whether it is lower than e^ t or not. 

Things do not get better if we set n + ^ 0. The plot in fig. I3.1bl was obtained using the 
same parameters as in fig. I3.1a[ but n + was set to 10% of n. In accordance to what we found 
in section T2.6I the lower order energy density rises; its minimum now lies above ej^ t . Also 
the higher order energy density increases and the fact that we cannot identify a minimum 
becomes even more distinctive. When going to lower values of n, the situation becomes even 
worse. 

Let us investigate the chiral limit. In section [231 we found that for = ao is fixed to 
7r/2. If this is still the case here we might get a result as the expansion is useful for large values 
of a . The plot in fig. ET2l shows e to t at n + /n = 0.4995 and e^ t , both at n 1 / 3 = 47.36 MeV 
and M„ = 0. 
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Figure 3.2: e tot and e£ s t at n 1 / 3 = 47.36 MeV, n + /n = 0.4995 and M w = 

Indeed, this plot suggests that the minimum of the energy density is located at ao = n/2. 
The mentioned value for n was found by the observation that, at ao = vr/2, the minimum of 
etot as a function of n + is located at n + = n/2. Then, the root of the difference e™ t — etot 
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at n + = n/2 and Qo = vr/2 as a function of n was found numericalljo to be 47.36 3 MeV 3 . 
However, at this moment we cannot make any statements about the behavior of etot f° r small 
values of ao an d we cannot exclude the possibility that the true minimum lies e.g. at ao = 0. 
We will make up for that later in section 13.1.21 where we will find that there is no minimum 
at small values of ao- Accordingly, we have found the critical particle density in the chiral 
limit to be 

n c = 106.2 GeV 3 or h x J z = 47.36 MeV. (3.35) 

However, as we noticed above, in this expansion we are not able to find the critical 
densities for = 40 MeV or = 135 MeV, since the minimum of etot then seems to be 
located in a region where our expansion fails. We therefore have to try to do the expansion of 
the fermionic energy ()2.34p explicitly for small values of ao an d first try it for ao ~ ^i 171 ]^)- 



3.1.2 Expansion for a ~ 0{ m N l ) 

Let us introduce a parameter rj of dimension of energy to write 

ao = . (3.36) 

m N 

Varying r\ in an appropriate range of 0(1) brings down ao to 0{rn N 1 ). (|3.36p implies 

P = a— + 0(m N 3 ) (3.37) 

WIN 

S = ^ + (m N 2 ) (3.38) 



and, when using (|3.2p and (j3.3[) . the internal square root of (|2.34p therefore may be expanded 
as follows: 



4(P ■ 5? + 9\{m N f$ 2 + (P ■ P) 2 } = $JgW +Pl + 0{ 



m 



N 



" 7 -,2^2 . „2 , /-or™- 2 



-^g 2 A V 2 +Pi + 0(m-'), (3.39) 

where in the last step we have expanded the second square root under the assumption that 
the terms 0(m N ) in the root are small. These are of the form 

2 2 4 

P3V V 



^N^A^+Piy mliiaW +P3)' 
Since p 2 <C m 2 N (see appendix [B]) the first term is indeed small as 

" 2 <i. 



9 A T + Pi 



The second term is larger than 1 only as soon as t] > QA m N •, if P3 = 0. For p 2 > r\ may 
even be larger. According to (|3.36p this would mean ao > qa = 1-26; we are however not 
interested in such high values for ao- Hence, in our regime the expansion of the internal 
square root should cause no problems. 



""^As it seems complicated to take the limits \xm n ,^ n /2 etot and lim ao -5.^/2 £ tot respectively we set n+ = 
0.4995n and ao = 7r/2 — 0.05 in the numerical evaluation 
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Plugging (|3.39p into (|2.34[) and replacing j3 and 5 now yields 



h± \ m% • • ^- : a^/g 2 A i] 2 +p 2 + 0(m N 2 ) 



mN \l 1 + iS" + T^T ± ^T\/9W+pI + 0(m- N 4 ) 



m N Am N m N 

p 2 a 2 



mN ( 1 + tt-t + ± 



2m^ 8m^ 27?!^ 
p 2 a 2 , a 



V^ 2 +p! + o(«^ 4 )) 



= ™ + 2^ + 8^ ± + *>" + <"*~ ). (3.40) 

The expansion of the square root in the second step is of course only allowed if 



m 2 N 



< 1. (3.41) 



As the dispersion relation ()3.40p now contains ps under a square root, the surfaces of 
constant energy in momentum space and hence the Fermi volume are more complicated 
than before. Accordingly, the integration over the Fermi volume in the calculation of the 
particle number N± and the total fermionic energy E^ ot will be a bit trickier. The calculation 
is performed in appendix O There, it turns out that for the calculation of n_ and e^ ot 
respectively we actually have to distinguish between two cases depending on the value of the 
spiral parameters a and rj, namely: 

• Case 1: a < IgAf] 

• Case 2: a > 2g^r]. 

The particle number density n_ and the energy density (the superscript (1) standing 

for case 1) are given by (jC.llj) and (|C.12[) respectively, while n_ and e[ { are given by 
()C.16p and ()C,17p , For the calculation of n + and ef ot no distinction is necessary and the 
results are found in (|CT2l]> and (ITI22]) . 

Unfortunately, the expressions for n_, and n+, which are functions of the Fermi 

energy Ep, cannot be solved analytically for Ep in order to express , e[^ t and e t ' ot in 
terms of n_ and n+. Accordingly, we cannot minimize the total energy with respect to the 
parameter a analytically in order to get rid of it. We rather have to do the replacements and 
the minimization numerically in order to find the total energy density etot as a function of 
ao at given values of n and n+ (or n_). This is done in a Mathematica program according 
to the following scheme: 

1. Split up a fixed particle number density n into two parts n_ and n + , such that n = 
n_ + n + . 

2. Fix the value of ao and hence of r\. 

3. Set a to a particular value, e.g. a = 0. 

4. Depending on whether a < Ig^j] or a > 2gAi] numerically solve either (jC.lip or (|C.16p 
for Ep and plug this into either (|C.12p or ()C.17p to find e tot . 



5. Numerically solve (|C.2ip for Ep and plug this into (|C.22p to find e 



tof 
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6. Calculate e to t = ^tot + e tot + e p> where e p is given by (|2,28|) . 

7. Increase a by a certain amount. 

8. Repeat steps 4-7 up to some appropriate value of a. 

9. From the set of all values of e to t obtained in such a way, choose the minimal one to be 
the total energy density for the given angle ao- 

10. Increase ao by a certain amount. 

11. Repeat steps 3-10 up to a certain value of ao- 

As we target to a critical number density from where on a spiral configuration will be favored 
we still need to identify the non-spiral energy density. A non-spiral configuration is given 
for a = 0. The plots in fig. 13.31 show e tot at fixed values of n and ao and at two different 
values for n_. They were obtained by accomplishing the above steps 1-8 and show what we 
suppose: The non-spiral configuration with a = is given for n_ = n+ = n/2. 




100 200 300 400 100 200 300 



a [MeV] a [MeV] 

(a) n_/n = 0.5 (b)n_/n = 0.65 

Figure 3.3: etot as a function of a at n 1 / 3 = 60 MeV, rj = 15 MeV and a physical pion mass. 



Hence, let us calculate e to t for o = and n_ = n + = n/2. Since a = < 2^7/ we have to 
find e^ t in the framework of the first case. For a = (jC.llj) reduces to 



(l) m N 
rvj = — 7T 

7T Z 



F 

J dE yj2m N (E - m N ) = ^(2m N (E~ - m N )fl 2 . (3.42) 



m N 



This may be solved for E F , and plugged into (|C.12j) . leading to 

e tot = m N n. + JL- (9vr V ) 1/3 . (3.43) 

For a = ()C,2ip reduces to exactly the same expression as (|3.42|) and (jC22|) to (|C.12p . 
except that we have to replace Ep by Ep and n_ by n + . e^ ot is therefore given also by (|3.43p 
with n_ — > n + . Finally, setting n_ = n + = n/2 the total energy density reads 

etot = e t ~ot + 4ot + e P = m Nn + j^- ( -vr 4 r?, 5 J - F 2 M% cos a , (3.44) 

which is minimized for ao = 0, such that we get exactly the same non-spiral energy density 
as in (|2.6ip or (|3.34p respectively. 
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Now that we have everything together we can ask for a critical particle number density. 
Based on the insight that at first order for physical pion masses n + will go to zero (see section 
I2.6p we now make the same assumption here and just take n_ into account. Later we will 
verify this assumption and in fact learn that it is not exactly true anymore. 

The "appropriate value" of a in above step 8 depends on the values of n and 77. It is 
clear that it must not be too large, otherwise the condition (|3.4ip will not be satisfied and 
the expansion (|3,40p is not valid anymore. Going to too small values of rj (for a given n) 
will manifest in the fact that etot will not show a distinct minimum as a function of a. This 
can be seen in the plots of fig 13.41 which where also generated numerically through the above 
steps 1-8. 




100 200 300 400 500 600 700 100 200 300 400 500 



a [MeV] a [MeV] 

(a) r\ = 30 MeV (a = 0.032) (b) r\ = 5 MeV (a = 0.005) 

Figure 3.4: e to t as a function of a at a density of n 1 / 3 = 60 MeV 



It is obvious that in the first plot (ao = 0.032) a clear minimum of etot around a = 200 MeV 
may be found. For values of a larger that 500 MeV etot drops again and in fact starts to diverge 
towards — oo. For too large values of a condition (|3.41|) is not satisfied and the expansion 
not valid anymore. In the second plot (ao = 0.005) etot starts to diverge for smaller values 
of a. In that particular case the drop starts so early that we are not able to identify a local 
minimum. We should therefore hope that the minimum of etot as a function of Qo around 
the critical density n c does not lie at too small values of ao, otherwise it can not be found in 
the framework of this expansion and we should do it for even smaller angles. 

Fig. 13.51 shows the same plot as in fig. 13. la|, i.e. the the lower and higher order energy 
densities as well as the non-spiral energy density at n 1 / 3 = 102.27 MeV, supplemented with 
the new energy density obtained through above steps 1-11. Unlike the previous higher order 
energy density the latter one shows an obvious minimum which, as it lies below the non-spiral 
energy density, allows us to state that we are already above the critical number density n c . 
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Figure 3.5: 1: First order energy density, 2: Higher order energy density of (|3.29j) and (|3,30p . 
3: Higher order energy density in the ao ~ ©(m^ 1 ) expansion, 4: Non-spiral energy density 
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In order to find n c we repeat the above 11 steps for a certain range of n, each time picking 
out the minimal value of etot (as a function of ao) and calculating the difference d between 
e^ t and the minimum (cf. the function d max (n) in section I2T0]) . The result is the plot in fig. 
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Figure 3.6: Difference between e^ t and the minimum of e to t as a function of n 

The root and thus the (alleged) critical number density is found to be located at 

n = 425.3 GeV 3 or n 1 / 3 = 75.20 MeV. (3.45) 

We should now check whether the assumption that also here for physical pion masses 
only EL contributes to the fermionic energy density is true. To do so we apply the above 
11 steps to the particle density n 1 / 3 = 75.20 MeV we just found but do not restrict ourselves 
to n + = 0. The plot in fig. 13.71 shows etot at n 1 / 3 = 75.20 MeV and n + = (continuous 
graph) as well as e to t at the same particle density but with n + /n = 0.001 (dashed graph). 
The horizontal line corresponds to e^ t . 
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Figure 3.7: The minimum of etot drops as n+ brings a contribution. 

In contrast to our assumption the energy density and in particular its minimum does not 
increase but drops for n + ^ 0. This differs from the properties observed earlier and shows 
that we actually have not yet found the critical particle number density. Since of course 

1/3 

ej^t is n °t influenced by n+ we have n c < 75.20 MeV. The problem is now the following: 
The lower n is chosen the more will move the minimum of e to t towards smaller values of 
ao- It then becomes numerically more and more difficult to identify a minimum. The lowest 
approximate value of n where one is able to find a minimum of etot which lies below e^t lies 
in the region around n = 398.9 GeV. For densities lower than this value the minimum seems 
to have shifted to ao = and to have approached e^. The plot in fig. I3.8al shows e to t at 
n = 399.1 GeV 3 and for n + /n = 0.07. We can identify a minimum slightly below e{^ t which 
shows that we are just above the critical particle density. On the other hand, the plot in 
fig. I3.8bl was created using n = 398.8 GeV 3 . At this density it is not possible to choose n+ 
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(a) n = 399.1 GeV 3 , n+/n = 0.07 (b) n = 398.8 GeV 3 , n+/n = 0.3 



Figure 3.8: e to t and ef^ t slightly above and below the critical particle density 

such that the minimum of etot lies below the non-spiral energy density; we are thus below 
the critical particle density. From this we conclude that n c is located somewhere between 
398.8 GeV 3 and 398.9 GeV 3 . We take the larger value as an upper bound: 

n c = 398.9 GeV 3 or n l J ?J = 73.61 MeV. (3.46) 

Compared to the first order result ()2.63p this is a 28% decrease of the critical particle number 
densityH 

In the same way as we found n c , we also determine n c (i.e. using M n = 40 MeV): 

n c = 172.1 GeV 3 or n l J z = 55.62 MeV. (3.47) 

In section [3. 1.1 1 we found the critical particle density in the chiral limit under the assump- 
tion that the minimum of etot is at ao = tt/2 and not in the region of small values of otQ. 
With the expansion for qo ~ 0(mj^) we can now close the gap and investigate the behavior 
of etot for small values of ao- The plot in fig. !3.9l shows etot as a function of Qo using the same 
parameters as in the discussion of the chiral limit in section f3 . 1 . 11 i.e. n 1//3 = 47.36 MeV, 
n+ = 0.4995n and M w = 0. 




0.05 0.1 0.15 0.2 

Figure 3.9: e tot at n 1 / 3 = 47.36 MeV, n+/n = 0.4995 and M n = 

From this plot we see that there is indeed no minimum of etot at small values of ao- Hence, 
we subsequently justify the result for n c found in section [3.1.11 

Apart from the chiral limit where ao = 7r /2, we found that, at the critical density, the 
minimum of etot seems to lie at ao = 0. This is a good sign. Whenever we figured out the 
non-spiral energy density e^ t we found that ao will go to zero. It is therefore peculiar that in 
the lowest order calculations (see section |2T6|) ao seemed to adjust itself to a non-zero value at 
the critical density. The transition from the spiral the the non-spiral configuration would then 
show an unnatural discontinuity. Here, on the other hand, the energy density "smoothly" 

2 As we want to make a comparison on the energy scale we compare n 1//3 rather than n. 
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approaches the non-spiral energy density when shifting towards n c from above. As a non- 
spiral configuration (so far) always came with an equal filling of the states (n+ = = n/2) 
we assume this to happen between 398.9 GeV 3 and 398.8 GeV 3 (for the physical pion mass). 

3.1.3 Expansion for a ~ 0(m~^) 

If this is true and the minimum of etot indeed approaches ao = when n \ n c , then one 
should verify whether the spiral parameter a goes to zero for ao -4 0. The plot in fig. I3.8bl 
gives a hint that this is in fact true since there the minimum touches ef^ and this happens 
at «o = 0. Additionally, the plot in fig. I3.1UI shows the parameter a minimizing etot of the 
first order calculation (see (|2.42p ) as a function of ao- It shows the described behavior (even 
though at first order the minimum of etot does not approach ao = at the critical density) . 
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Figure 3.10: a m j n at first order with ra 1 / 3 
70MeV and n + = 
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Figure 3.11: a m j n at next order in the ao ~ 
0(1) expansion with n 1//3 = 70MeV and 
n, = 



On the other hand, the plot in fig. l3.11l shows the minimizing a at next order in the ao ~ C(l) 
expansion, i.e. a + of (|3.3U|) . Here, a — > oo as ao — > 0, which would mean that the spiral 
would rotate faster and faster as we go to smaller values of ao- However, this would be bad 
as the energy expansion for too large values of a fails, or, equivalently formulated, as a too 
quickly rotating spiral would imply too large spatial derivatives and therefore too large pion 
momenta for the ChPT to be valid. 

In order to verify our supposition that the spiral disappears for ao — > 0, let us therefore 
expand the fermionic energy at even smaller values of ao than we did it previously, namely 



for ao ~ 0(m N ). The expansion of ()2.34p with ao 



K 2 /mjf 



(analogon to (|3.36p ) reads 



E± = m N + h ± 



2mjv 8mjv 2niN 
which corresponds to (|3.4U|) with rj = 0. As 



+ 0(m- 3 ), 



(3.48) 



. a k 
o = — cos 
2 m 



2 --- + 0(m^ 4 ), 

N Z 



we replace a/2 by 5 in (|3.48|) : 



E±=m N + h ± 

2m7v zttt-tv tum 



(3.49) 



We now first treat E + . We note that surfaces of constant energy in momentum space are 
symmetric with respect to and that we hence can integrate over p% > only and then 
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multiply by two. Thus, we can drop the absolute value in (|3.49|) to get 

E+ = m N + -L- (pf + p\ + (p 3 + 5f) . (3.50) 

zmjv 

From this we see that a surface of constant energy is a sphere centered at P3 = —5, however 
only for p% > 0. For p% < we have the mirror image. Fig. 13.121 shows a cut of a surface of 
constant energy in the pi — ps plane. 




Figure 3.12: Cut of a surface of constant energy for E + in the p2 — P3 plane 



Shifting the integration variable P3 to p' 3 = p% + 5 allows us to integrate over the upper part 
of the sphere using ordinary spherical coordinates, where the polar angle 6 runs from to 
^ = arccos bjp. The absolute value p of the momentum runs from 5 to some momentum pf+, 
which is the momentum reaching from p% = —5 to the surface of the Fermi volumeU These 
considerations are only true if pf+ > 5 (otherwise we will get no contribution from E + ). For 
large particle densities and if we restrict ourselves to small values of a this is satisfied. For 
the particle number we then get (note the additional factor of 2) 



iV 4 



<0 



PF+ 2tt £ 

dpp 2 J dip J d6 sin( 




PF+ 



dpp 



dz 



PF4 



^2 / ''I'l'' 



I- 6 - 

P 



3tt 2 



35 
2p F ^ 



+ 



5 

S 3 
2p 3 



S/p 



(3.51) 



where in the first step we performed the variable transformation z = cos 6. Analogously we 
find the total energy due to E + states: 



2tt 



mNN+ + 10^7 ( 4+ " I 



o o 

.5 



dpp 2 dip d9 sin# ( + 



1 



2m 7v 



V 



(3.52) 



Next we consider Due to symmetry we may take again only p% > into account. 
The dispersion relation may then be written as 



1 



m N + {pi +p 2 2 + (p 3 - 5f) 



(3.53) 



3 Since the Fermi volume is not a sphere we cannot call this momentum Fermi momentum. Nevertheless 
we denote it by pf- 
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which shows that a surface of constant energy is again a sphere but this time centered around 
P3 = 5. For £>3 < we have once again the mirror image. Fig. 13. 131 shows again a cut of such 
a surface in the p2 — P3 plane. 



s 









Figure 3.13: Cut of a surface of constant energy for £7_ in the p2 — P3 plane 



After shifting the integration variable P3 top' 3 = P3 — S we again use spherical coordinates. Let 
us denote the absolute value of the momentum that reaches from p 3 = 5 to the Fermi surface 
by pf~- Again under the assumption that pf- > <5 we now have to split the integration into 
two parts: A first part where p < 5 and where the polar angle 9 £ [0, ir] and a second part 
where p > 5 and where 9 runs from to £ = ir — arccos 5/p: 

3 S 2tt it Pf- 2tt C 

iV_ = 4 (J^j J dp p 2 J dtp J d9 sm9+ J dpp 2 J dp J d9 sin# 
5 

L 3 3 L 3 1 

= 16)r (^) T +ftr (£) I dpp2 1 dz 

S -5/p 
Pf- , N 

2L 3 5 3 L 3 f 2 / 5\ L 3 p 3 F _ 35 5 3 \ , nrA . 

+ ^ P ! + - 1 + — - -r . (3.54) 



37T 2 7T 2 J \ p J 3ir 2 I 2p.p_ 2p|,_ y 
Similarly the total energy due to 2?_ states is calculated, the result being 

E « - mwJV - + loJ^ + 1 { " 4 "- - t) ■ < 3 ' 55 ' 

Dividing (|3.5ip , (|3.52p , (|3.54p and f|3.55j) by L 3 we may write the particle number and energy 
densities in summary: 

e f ot = m N n ± + U ± =F 5 -5p F± ± ^) . (3.57) 
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(|3.56|) may be solved for pf±] the expression reads 

PF± = \ [k{n±,T5) ±5}+ S —— , (3.58) 

where 

k(n, 5) = (l2mr 2 + <5 3 + 2iry/36n 2 ir 2 + 6n<5 3 ) . (3.59) 

Plugging (|3.58p into (|3.57|) and substituting ra_ = n — n+ finally allows us to express the total 
fermionic energy density ej = ef ot + e^ ot in terms of of n and n + . The plot in fig. 13. 141 shows 
the total energy densit)|j e to t = e/ + e p at n 1//3 = 73.61 MeV as a function of n + and 5. 




Figure 3.14: The total energy density at very small values of ao and n 1//3 = 73.61 MeV 

The plot shows a similar behavior for different values of n and confirms what we wanted to 
show: For very small values of ao the parameters 5 and hence a go to zero and the spiral 
disappears. 



3.2 Higher Orders in the Baryonic Sector 

Having gone away from the chiral limit and therefore including chiral breaking terms in the 
pion sector already in section 12.61 we should also consider these in the baryonic sector. This 
leads to the next order of the fermionic Lagrangian, which is given by (|1.65|) . In our spiral 
configuration we have 

[uo,^] = (since uq = 0) 

[Ui,Uj] = [Ci,Cj] = Ci, a T a ,Cj, b T b = [ci^T 2 ,Cj^T 2 ] = fiifij [t 2 , T 2 ] = 0, 

i.e. the term proportional to C4 disappears. Furthermore, if we switch off all external fields 
except s = M. = diag(m u , nid), X of (|1.54p becomes \ = The second order baryonic 

Lagrangian then reads 



C {2) N = 2B Cl {u ] Mu ] + uMu)^^ - -%-(u iUl ) (WD% + h.c.) + — Imu*)^. (3.60) 

4m^r ' ' 

Let us introduce the following notation: 



t 1 = 2Bci{u ] Mu ] + uMu), tl = %-iuiUj), t 3 = —(u iU l ). (3.61) 

J 4m^ 2 



4 At this order e p = -F 2 Ml + 0(m~ r 
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Using (|1.60|) , (|2.2|) and the fact that the connection = cL is antihermitian (see 12. 8p we 
find 



WZ>% + h.c. = ^D l D j $> + (tfZWtf) t = vpZ^D-% + Nft'y D i D i & 

= * (iW + D^iA ^ = q, + d i)(Qj + d i) + (#/ _ _ rfi)] ^ 

= ^ (2d i d j + d l d? + d j d l ) ^ = ^ (2d i d j + {d\ d j }) ^. 
Remembering (|2.30p we have d % = —i5 l T 3 and therefore 

{d\d j } = -5 i 5 j {r 3 ,r 3 } = (3.62) 

where of course the 2x2 unit matrix is not written (exactly as with the term 2<9 J <9 J ). We 
can now write the Lagrangian as 



t 



C^ N = (t 1 + t 3 ) tftf + 2tj^ (d l d j - 5 l 5 j ) 



(3.63) 



We now go on as in section I2.1.1| i.e. extract the single-particle Hamiltonian via Euler- 
Lagrange equations. Its individual parts read (using the full Lagrangian Cf = C^ N + c!^ N ): 

•(1) flr( 2 ) 



dC 



f 



=0, see 

z. fern 











dC f dC 



(i) 

ttN 



+ 



dC 



(2) 
:N 



f L +(t 1 + t 3 ) ^ + 2t% - 5 l 5 j ) f . 



d^S> 

The term dC^ N / d^f has already been calculated in section 12.1.11 namely in (|2.5p . In this 
term the sum ^^d^ appears. This sum was split up into the temporal and spatial part, the 
spatial part was brought to the other side of the equation of motion and the whole equation 

was multiplied by 7 in order to have the combination id^ isolated on one side of the 

(2) 

equation. Since there is no additional time derivative originating from C w ^, these steps can 
immediately be repeated: 

EOM: = i^d^ + if d^ - mjv* + y^c^ + (i 1 + t 3 ) V + 2t?- (d i d j - <f (P) f 



9A 
2 

9A 



%1% - i-y% + m N - ^ 7 SsCi - {t 1 + t 3 ) - 2t% {&& - 

fj 5 c l -{t 1 +t 3 )-2t* 1 {dW -5 l V) 



H^ - (t 1 + t 3 ) 7°^ - 2t\f - &&) ^. 



H\ is the first order Hamiltonian of (|2.7p . while the rest can now be identified as H%. Using 
= d i (-ip'^z) = -ipW^ = -pVf, we get 



Ho 



(t 1 + t 3 ) 7 + 247 (pV + 5 l V) 



(3.64) 



Comparing the 7 matrix structure of H2 with Hi we see that H2 contains only terms propor- 
tional to 7 (actually diag(7°, 7 )). If H\ = 017° + gT is an operator (with T being a linear 
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combination of 7 matrices or products thereof, g denoting the set of coefficients) leading to 
some eigenvalues f{a%,g), the eigenvalues of a new operator H = H± + 027° are easily found: 

H = a^ + gT + a 2 j° = (ai + a 2 )7° + 5 r — > f(a 1 + a 2 ,g). (3.65) 

In Hi it is the mass term being proportional to 7 . Hence, in order to find the new eigenvalues 
of the full Hamiltonian H = H\ + H 2 we only have to replace 

m N — ► m N + 2t% (pV + S l 5 j ) - (t 1 + t 3 ) (3.66) 

in the old eigenvalues of H\ . 

Let us calculate the remaining quantities: As Ai is diagonal it commutes with u and 
and we can write 

t 1 = 2B Cl (MU j + MU). 

This trace we have already calculated in (|2.26|> . Assuming again the isospin limit m u = = 
m we therefore get 

t 1 = 8Bcim cos «o = 4M 2 ci cos ao, (3.67) 

where in the last step we used (jl.42p . 

The trace {uiUj) we calculate using the constant vielbein of (|2.23p . (|2.8p and the definition 
(I239|) . For tfj we then get 



," C 2 / \ C 2 / 2 2\ C 2 Pi/3j C 2 Pi(3j ,o,<o. 

= -4^ te) = "4r< (C *< 2r C - 2T } = = (3 - 68) 



In the replacement (|3.66p we now can write 



N 2 



^=-2^-()9-*) • ( 3 - 7 °) 



4-2 xid — C ?_ 

Finally, we calculate the trace {uiU 1 ) and find 



t 3 = y(ci^) = y& i - r, ,V - -r, i, >r . (3.71 ) 



The replacement ()3.66p now reads 



C2 

— ► mjv 5- 



+ c 3 p 2 - 4M 2 a cos a . (3.72) 



After having performed this replacement in the energy expression of (|2.34p we expand it 
in powers of l/rrjjv up to order 1/mjv which yield^l 

V 2 S 2 1 
E± = m N + + + c 3 /3 2 - 4 Cl M 2 cos a ± - 5A /3. (3.73) 

Since at this order the energy only depends on p 2 we again talk about a Fermi sphere in 
momentum space and accordingly can make use of the formula (|2,37p and just slightly modify 
(12~39D to 



efot = [m N + + c 3 f3 z - 4 Cl M^ cos a ± -g A fi n± + V , — . (3.74) 



(37r 2 n±) 5 / 3 



5 i.e. we again consider the lowest order. Even though in this expansion we found the peculiar behavior 
that the minimum of etot does not approach ao = as n \ n c , we still can compare the results with (|2.55p 
and (|2.63|) in order to get an impression of the effects the higher order terms in the Lagrangian have. 
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Correspondingly, the total fermionic energy density now reads 

e / = e tot + e tot 

( , 52 l_ R2 a \ , , * , (37T 2 ) 5 / 3 / 5/3 , 5/3 

= m JV + 7; V C3P - ■iciM^ cos a n + — - (n + - n_) + — — ^ n / + n_' 

\ 2mN J 1 LUTT TYl]\[ \ 

(3.75) 

Since the low energy constants Cj are of order m^ 1 (see section |1.8.3|) the expression (|3,75p 
shows us that, in terms of an expansion in powers of , we are actually not examining 
higher orders here, but rather completing the first order expansion. 

In the pionic sector we stay at order 0(p 2 ) and therefore again use (|2.28|) for the pionic 
energy density. After having replaced f3 and 5 in terms of a and ao and setting n_ = n — n+ 
the total energy density is 

/ a 2 cos 2 ao 22 * r i \ 9 A 
etot = rriN H h £30 sm ao — ^c\M cosqo n asmao n — 2n + ) 

\ &TYI ' / 2 

+ \2m N ( n + 3 + n+)5/3 ) + T ^ Sin2 °° " COS a °' (3 ' 76) 



The value of a that minimizes etot is found to be 

2rri]ygA sincere — 2n H 



n cos 2 «o + 4mAr(F 2 + 2c3n) sin ao ' 



(3.77) 



which, after plugging into (|3.76p leads us to the total energy density as a function of n, n + 
and ao- 

3.2.1 Chiral Limit 

Even though we basically included higher order terms of the fermionic Lagrangian in order 

(2) 

to take into account the effects of chiral symmetry breaking, C L still has terms that do not 
explicitly break chiral symmetry. Let us investigate their influence by considering the chiral 
limit, such that we can compare the results with those obtained in section [231 

Varying e to t with a replaced by a m ; n we find in complete analogy to section 12.51 the 
condition sin2ao = 0. To decide which one of the two solutions, ao = or ao = vr/2, 
minimizes the energy density, we calculate the second derivative of etot at a o = 7r/2 and get 
an expression very similar to that of (|2.45j) : 



<9 2 e to t 



g\n{n - 2n H 



«o=f 16m 7V (F 2 + 2c 3 n) 2 - (3 ' 78) 



Again this is positive for < n + < n/2 and we conclude that, as before, in the chiral limit 
a o = 7r /2, i.e. the pion vector moves in the equatorial plane. The expression for the total 
energy density then simplifies to 

e » - m " n - 8( f/+L 3 „) '" - 2n +> 2 + IIS ("+ 3 + <» - n +> 5/3 ) ■ < 3 ' 79 » 

We vary this expression with respect to n + , 

aetot _ (3vr 2 ) 5 / 3 / 2/3 _ , )2/3 \ _ 9 2 A (2 _ n) l Q f3 80) 

dn+~&r*m N \ + [ +) > 2(F 2 + 2c 3 n) ^ + U) ~ U ' (6 ' m) 
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and find (cf. section |2~5"]) that for small n this equation has the only solution at n + = n/2, 
which, because of (|3.77p . means a non-spiral configuration. Increasing n above the critical 
particle density, the second derivative of etot with respect to n + at n + = n/2, 



<9 2 etot 



dn 2 + 



2(3vr 2 ) 5 / 3 /2\ 1/3 g\ 



97r 2 mAr \n J F 2 + 2c3 



n 



(3.81) 



gets negative and the extremum at n+ = n/2 becomes a maximum i.e. the minimum lies 
somewhere at n + ^ n/2 and the spiral configuration will be favored. The critical particle 
number density n c is the density where the second derivative becomes zero. For C3 = — 3Am7} 
(s. section 11.8. 3p this is the case at 



79.2 GeV 3 or nj /3 = 42.95 MeV (3.82) 
Using C3 = — 4.2m^ r 1 the critical density even drops more: 

75.5 GeV 3 or nJ /3 = 42.26 MeV (3.83) 



n r 



(2) 

Hence, due to the additional terms of the higher order baryonic Lagrangian C ^, compared 
with (|2.55j) the critical particle number density in the chiral limit drops by about 8 — 10%, 
depending on the exact value of C3. 



3.2.2 Back toM^O 

Let us now go back to (|3.76p . Since ao 7^ vr/2 for M n / we have to deal with both n+ and 
ao- We treat this problem as we did in section [2T6l We notice that due to (|3.77p a non-spiral 
configuration is given for n + = n/2. This leads to 

(37r 2 n) 5 / 3 

e tot (n+ = n/2,a ) = m N n + \ ' M 2 (F 2 + 4cm) cos a , (3.84) 

2°' ,i OTr z m]sr 

which is minimized for ao = (if c\ = — O^mJ^ 1 (— 0.6mJ^) the expression F 2 + 4cin > for 
n 1 / 3 < 121.9 MeV (139.6 MeV)). 

As in (|2.62|) we now define the function d(n + ,ao) as the difference between the energy 
densities of the non-spiral and the spiral configuration. Plotting then d(n + ,ao) for different 
pion masses M n we find a behavior very similar to that found in the lower order calculation 
in section 12.61 With increasing pion mass the critical parameter n c + (i.e. the value of n + for 
which d{n c + , ao) = for exactly one ao) goes to towards zero (cf. figure [273]) . For pion masses 
M w > 20 MeV we can restrict our analysis to the parameter subspace n + = 0. We then 
search the maximum values of d(0,ao) for a given set of values for n, defining the function 
d max (n). The critical particle number density n c we are looking for satisfies d m&x {n c ) = 0, 
because then the energy density of the spiral configuration is equal to that of the non-spiral 
configuration. The plot in fig. 13.151 shows the function d max (n) for both c\ = — 0.9m^ , 
C3 = — 4.2m^ r 1 (continuous graph) and c\ = — 0.6m^ , C3 = — 3.4m^ (dashed graph) at the 
physical pion mass. The two sets of low energy constants Cj are the ones given in section 
11.8.31 The critical densities for c\ = — 0.9m^ , C3 = — 4.2m^ r 1 are numerically found to be at: 



n c 



410.8 GeV 3 or n^ 3 = 74.34 MeV (3.85) 
196.5 GeV 3 or ¥lJ 3 = 58.14 MeV, (3.86) 
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Figure 3.15: d maiX (n) with the physical pion mass and for c\ = — 0.9rrt^ , C3 = —4.2m 
(continuous) and c\ = — O.Gm^ 1 , C3 = — 3.4m^f (dashed) respectively 



while for c\ = —0.6m N , C3 = —3.4771^ we find 

n c = 465.9 GeV 3 or n^ 3 = 77.52 MeV (3.87) 
n c = 213.4 GeV 3 or n l J 3 = 59.75 MeV. (3.88) 

Compared to the critical particle number density we found at lower order in section I2.6( for 
M w = 135 MeV this means a 24 - 27% decrease and for M w = 40 MeV a 15 - 17% decrease 
of the numerical values, depending on the values of the low energy constants c\ and C3. 



Chapter 4 

Conclusion and Outlook 



In this thesis we have shown that in the framework of ChPT the pionic vacuum may change 
from a homogeneous to a spiral phase if the density of fermions is chosen to lie in the region 
of nuclear matter density. The numerical values for the critical fermionic density are collected 
here: 

— 1 (2) 

In the expansion of the fermionic energy up to 0(m N ) but without the terms of C^L 
(see sections 12.51 and I2.6P we founcQ 

nj/ 3 = 46.81 MeV n l J 3 = 70.16 MeV n]J 3 = 102.27 MeV. (4.1) 

(2) —1 

Including also the chiral breaking terms of C^jy bringing additional contributions of 0(m N ) 
(see sections 13.2.11 and I3.2.2P , the values change to 

h l J 3 = 42.95 MeV n l J 3 = 59.75 MeV n]J 3 = 77.52 MeV, (4.2) 

when using c\ = — O^m^ 1 , C3 = — 3Am]^. Finally, going to higher orders in the expansion 
of the fermionic energy (see sections 13.1.11 and 13. 1.2[) led us to 

h l J 3 = 47.36 MeV n l J 3 = 55.62 MeV n]J 3 = 73.61 MeV. (4.3) 

The influence of the uncertainties in the low energy constants was found to be relatively small 
(see section I2.5P . 

We have several possibilities of comparing the results originating from different order 

calculations. Let us first take a look at the result obtained using the physical pion mass. 

(2) 

When excluding the terms of Aliy and going from the first order expansion of the fermionic 

energy to next order, the critical density drops from 102.27 MeV to 73.61 MeV. This is a 

28% decrease, which is quite large. Staying at the lowest order of the fermionic energy 

(2) 

expansion but including the terms of £1^ the critical density drops to 77.52 MeV, which is 
a 24% decrease. However, as mentioned in section 13. 2| including the terms of the next order 
baryonic Lagrangian was in principle completing the expansion of 0{rn^). It would therefore 

also be interesting to examine both steps at once, i.e. go to 0{m~^) in the expansion of the 

(2) 

energy eigenvalue and to include the terms of C n ^. 

Let us have a look at the results that were obtained using M n = 40 MeV. The two ways of 
going to higher orders just mentioned imply corrections of 21% and 15% respectively. These 
are more moderate than in the case of the physical pion mass. Finally, in the chiral limit 

going to next order in the expansion of the fermionic energy yields an increase of the critical 

(2) 

density of 1% while including C n ^ causes a decrease of 8%. 



1 As a reminder: The symbol o denotes the chiral limit while a bar indicates a pion mass of 40 MeV. 
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Hence, the smaller the pion mass the smaller are the corrections. The corrections in the 
case of the physical pion mass should remind us to stay skeptical and that the numerical 
estimates we make are maybe not very reliable. However, they are nevertheless not too large 
and we may get confidence that the qualitative picture is correct. 

The spiral configuration we investigated in this thesis was an ansatz that allowed analytic 
solutions of the Dirac equation. It is however by no means clear whether there is not another 
inhomogeneous pion configuration which could lower even more the total energy density of 
the many-particle system or which could happen at even smaller particle densities. Hence, 
the numerical values we found in this thesis should be understood as upper bounds of critical 
densities for a transition to an inhomogeneous pion configuration to happen. 

What is the physical meaning of a spiral configuration we investigated? If we believe 
our calculations, then already in ordinary nuclear matter the pionic vacuum could differ 
from the vacuum outside the volume of dense matter. It is therefore not possible to apply an 
appropriate gauge transformation in order to transform the pion field everywhere to U(x) = 1 
and to establish the ground state of no pions. Hence, observed from outside, inside the 
dense matter pions could appear, which is known as pion condensation. Different from free 
pions, which decay in a rather short time, our configuration is static and the pions should be 
considered as a collective phenomenon rather than free particles. A comparison for this might 
be the neutron; in the free case it decays after approximately 15 minutes, but as soon as it is 
enclosed in an atomic nucleus it turns out to be stable. Pion condensation in dense matter 
has been studied by several authors using different approaches, see e.g. [251 E3 Ell E21 E3]- 

Although the structure of the interior of neutron star^l is not yet known in all details, 
there are models that include the phenomenon of pion condensation inside these superdense 
star remnants (see e.g. [23 E2 129 ESI G2B EH 133] ) ■ 

The results in this thesis might hence give some alternative ideas in an approach towards 
an explanation of pion condensation in nuclear matter and neutron stars. 

An interesting thing that could be done in further considerations is the analysis of higher 
order terms in the mesonic Lagrangian (see section II. 7p and going to even higher orders 
in the baryonic Lagrangian. So far, the calculations at 0(p 2 ) in the mesonic sector were 
only tree level calculations. This did not make it necessary to renormalize the low energy 
constants; we could take them as fixed and more or less precisely known values. Going to 
0(p 4 ) would however imply one- loop calculations with £2 and tree level calculations with £4. 
It is an interesting question how one could then calculate the energy of the system. Also if 
one goes to higher orders in the baryonic sector one will have to deal with loops. It is then 
by far not clear whether a Hamiltonian as we constructed it is still a well defined object. 
An alternative way of finding the energy of a nucleoli in the background of pions could be 
to consider its propagator. This could be calculated using a well defined Lagrangian and it 
should be possible to make a conclusion for the nucleon energy via the pole of the propagator. 

Another direction worth exploring is to switch on electromagnetic effects, i.e. to include 
virtual photons. Beside the quark mass matrix the quark charge matrix is an alternative way 
of breaking chiral symmetry explicitly and it would be interesting to find out the order of 
magnitude of the effects arising therefrom. 

In the case of antiferromagnets the transition from a homogeneous to a spiral phase of 
the staggered magnetization implies a spontaneous breaking of translation symmetry on the 
lattice. The consequence is a further Goldstone boson (helimagnon) . In the case of QCD 
we can expect a similar phenomenon. It might be interesting to study the properties of the 
resulting additional Goldstone bosons. 



2 For a short discussion on neutron stars and an estimation of their densities see appendix [Dl 
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Appendix A 

The Spiral Configuration 



In this appendix we derive the static spiral configuration ()2.1ip from the condition of having 
a constant connection and a constant mass term in the Lagrangian. For this purpose we use 
a different parametrization of the matrix U than in (|1.30p . namely 

1 



U(x) = — [<7(s)l+*T-7r(s)], a[x) = ^F 2 -tt 2 (x), (A.l) 
F 

where the three Hermitian fields 7r* describe pion fields. cr(x) is chosen such that UU^ = 1. 
For the sake of simplicity let us absorb F into the fields, such that we can write 



U{x) = a{x)t + m-V(x), a(x) = y/l - n 2 (x). (A.2) 
We require the mass term (Ai (U + U')) to be constant: 

U + = 2a(x) • 1 = const. => a = const. => tt 2 = const. (A.3) 

In order to calculate the connection Tj we need the matrix u, such that u 2 = U. We 
parametrize it as follows: 

u(x) = f ■ 1 + ipV. (A.4) 

TherQ 



f + gigi T i T j =f + g*g* = (f + g 2 ) t = 1; 



from which we find 



f 2 + g 2 = l => / = Vl-9 2 - (A.5) 

We calculate 

u 2 = f 2 + 2i/<?V - ffVrV = / 2 - 5 V + 2i/^V = (/ 2 - 5 2 ) 1 + 2i/^V, 

such that the requirement u 2 = U yields 

f 2 -g 2 = a (A.6) 
= tt*. (A.7) 

Let us now calculate the connection Tj (without external fields) according to (|1.6ip . The 
ingredients are easily found to be 

u^diu = fdif + ifr'ihft' - iu'T'ilf + ij kl T l <?d ig k + gid ig i 
udiv? = fdif - ifridigi + ig j r^d i f + i^ kl r l g^d t g k + <M</. 



1 using rV J = ie ljk T k + 8 lj t and, due to (anti-)symmetry, g t g j e' tjk = 0. 
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Using (|A.5p we have 

fd i f + g j d i g1 = ±d i (f 2 + g1gi) =0, 

1 ^ v ' 

1 

such that the connection may be written as 

T 4 = ie>'V%V. (A.8) 

Plugging ()A.5p into (|A.6|) yields 1 — 2g 2 = ct and since o = const, we conclude that g 2 = const, 
and therefore / = const. For this reason, inserting (|A.7j) into (|A.8j) leads to 

r* = i(? 2] di Tf = Tp^d^-r 1 . (A.9) 

If we now require the connection to be constant, we consequently must have 

€ jkl 7r j diTT k = const., or 7T x diir =: di = const. (A. 10) 

Due to (|A.3p we may write 



7T = I TTi, 7T 2 , a/1 - 7rf - 7r| 



diir = dim, diir 2 , -====- 



di = 7T X diTT = I - — (7Ti9j7ri + 7r 2 <9j7r 2 ) - 7r 3 5i7r2 , TT 3 diTTl + — (vri5j7ri + 7T2<9i7r 2 ) , 



The relation 7r 3 = y^l — nf — t: 2 allows us to simplify d\ and d 2 : 

d\ = — \\ 7r 2 7ri5i7Ti + 7T2«9jvr2 + 7^(9^2 ] = — df - ^—^ (A. 11) 

7T J \ v „ ' / 7T 3 7T 3 

(l-7rJ)9i7T2 

d* = — I tti^tti + vrf d^i -HrjTrato | = — df + — , (A. 12) 

^3 y „ ' J 7T3 7T3 

(1— 7r|)9i7ri 

where in the last step we inserted (if = it\d{'K2 — (jA.lip and (|A.12p are equivalent to 

dun = d 2 7T 3 - d^ 2 (A. 13) 

dnr 2 = dfiri - d}ir 3 , (A.14) 

which build, together with 

df = 7Ti3i7T2 - 7T2<9i7Ti, (A. 15) 

a system of differential equations for the 71*. Inserting ()A.13P and ()A.14p into ()A.15P yields 
df = df tti + d?7r| - d- 7rivr 3 - d, 2 7r 2 vr 3 =>- dj 7Ti7r 3 + d 2 7r 2 7r 3 + df (l - vr 2 - vr^) = 0. 
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The term in brackets is 7r|, such that the equation reads 

vr 3 (djiri + df 7T 2 + d-7r 3 ) = 0. 

This equation and thus (|A.15j) is trivially satisfied due to the relation 

??• 4 = 0, (A.16) 

which clearly holds since di = n x 9j7r. Thus, the two equations (|A.13p and (|A.14|) suffice 
to solve the problem. The di are 3 vectors in the 3-dimensional isospin space. They cannot 
all be linearly independent, otherwise relation (|A.16|) could not be satisfied. Hence, at least 
one of the dj must be linearly dependent from the others. However, if 7? is orthogonal to 2 
independent constant vectors and tt 2 = const., then also n is a constant vector. This would 
be trivial. Consequently, all 3 vectors di must be linearly dependent, i.e. they all have the 
same direction. We can choose the coordinate system in the isospin space such that they 
point into the third direction; di = (0,0, aj), where cij = \di\. Then, through relation ()A. 16|) . 
we immediately conclude that 

7T 3 = 0, (A.17) 

while flXjH) and (lATTil) simplify to 

di-Ki = -ai7r 2 (A. 18) 

9i7T2 = Oi7r 1 . (A. 19) 

Since aj is a constant, the solution of this system of differential equations reads 

tti(x) = Acos(aiXi + tpo) (A. 20) 

7r 2 (x) = Asin(aiXi + <po). (A. 21) 



From a = v/l — n 2 (x) it is clear that ir 2 (x) < 1. On the other hand we have 

7f 2 (x) = TT? (S) + 7rl(x) + 7r|(x) = 

showing that we must have — 1 < A < 1. This allows us to parametrize ^4 as 

yl = sina - (A. 22) 

According to (|A.2p the matrix f7 (x) thus finally reads 

J7(x) = cos aol + i sin ao (cos (pfflr 1 + sin (p(x)r 2 ) , (A. 23) 

where 

9?(x) = a ■ x + ip . (A.24) 

Hence, through the requirement that the connection and the mass term are constant we 
are led to the spiral configuration (|2.1ip . The vielbein can then be made constant by an 
appropriate gauge transformation, as we did it in section 12.2.21 



Appendix B 

Analysis of the Energy Expansion 



Here we want to go through the expansion of the fermionic energy (|2,34p in order to find 
conditions for the expansion to be true. Let us first define the vector 

n = \gA&. (B.i) 

With this, the energy eigenvalues of (|2.34|) may be written as 



E± = W m 2 N + p 2 + 5 2 + 7 2 ±2^m 2 Nl 2 + (p- 5) 2 + (p- j) 2 , (B.2) 
or, when using f|3.2j) and (|3.3|) . 



E± = \jm 2 N + p 2 + 5 2 + 1 2 ± 2^m 2 Nl 2 + p 2 (6 2 + 7 2 ). (B.3) 

Now we go through the expansion f)3 . 1 [) in more detail. First we process the internal square 
root: 



m 2 Nl 2 +p 2 (5 2 + 7 2 ) = rriNl Jl + ^2 + 7 2 ) = m Nl (l + - ) + 0(x 2 ) 



2 

™^v7 + + 7 2 ) + 0(m^), 



2?Tljv7 

where we assumed that 2 <C 1; 

pp 2 . A. ^ 1 



2 9' 2 

For chiral perturbation theory to be valid we must have p% <C mj^. p^, can at most be pf and 
at particle densities we talk about pp is indeed lower that mjv, according to the lower order 
result (|2.37p . Consequently, we require 

# _ pi * 2 ^ , . ^ 2 



«i -^<c»d). 



09 9 9 ' 9 rv^ 

mjyj z mjy 7 7 



Using (fBTTj) . (jZMD and QZHD this condition reads 



COS 2 «o 



2-2 ~ 

5^ sm z q 



<0(1). 
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For large values of ao this is always satisfied and for small ones cos Qo ~ 1 , such that we find 

5 isin 2 a >l (B.4) 

for the above expansion of the internal square root to be allowed. If this is satisfied we can 
then handle the outer square root in the same manner: 



E± = \ m 2 + p 2 + 5 2 +7 2 ±2mjv7 ±-^-(<5 2 + 7 2 ) 

V m Nl 

= m N Jl + J?+f T + J^±^± ^(S 2 + 7 2 ). 
y m z N m z N m z N m N 7171% 



Here, we have to take care of the fact that x contains a term ~ m ^ ■ We therefore have to 



expand up to x 3 producing terms ~ m N 2 and ~ m N 3 which yield terms ~ and 
after multiplying the prefactor tun back. So we use 



2 s 

JT+x~=l + -- — + — + 0(x i ), (B.5) 
2 8 16 v ; v ' 

which is again allowed for x<l only. For a not too large the only term in x that may cause 
problems is the last one, which blows up for small values of ao- So we must demand 

^.+4=4-— +4- JL <i =► — sew, 

7S m N m z N m N >y my m N m N ^ 

small small 

which is equivalent to 

~ — < I- (B.6) 

LrrtfqgA sin «o 

For large a also the other terms in x (except of the first one) may cause problems. Let us 
assume that we are allowed to do the expansion. Then, from x 2 we only pick the terms up 
~ 3 ) 

x 2 =±^±^±'4 + A 4 + o[ m - N % 

m% m% m% m A N 



to O(mJ), i.e. 



while from x 3 we have only one term ~ rn N 3 : 



x 3 = ±4r + °( m ~N 



m 



A' 



Using (|B.5p we are then led to 



IP 2 1^ 7 3 7 2 ± 7 3 
2m 3 N 2m 3 N 2m 3 N 2m 2 N 2m 3 N 



P 2 , 62 , , P3 /r2 , 2 W 7P 2 _,_ 7<^ 2 



m7V 2mjv 2?rijv ^ 27771^ ^ + 7 ) T 2 m 2^ ^ 277^ 
which is exactly the expression (|3.1|) . 



Appendix C 

Integration over the Fermi Volume 

in the ao ~ 0(m^r ) expansion 



Here, we perform the integration over the Fermi volume for the dispersion relation 



.2 „2 



which results in the expansion of the fermionic energy for ao ~ 0(m~^). 

Let us first find iV_. The trick is to extend the momentum integration from the Fermi 
volume to arbitrary momenta by introducing @(Ep — E-(p)), where Ep is the Fermi energy 
of the E- states. Next, we insert a 1 by using the property 



oo 

J dE5{E-E_(p)) = 1. (C.2) 



Hence, 



Fermi 
volume 



o OO 



2vr 



J dE J d 3 p 5{E - E_{p))Q{E F - £7_(p)). 



For E < the delta function will always be zero as E_{p) > 0. If E > Ep we should also 
have E-(p) > Ep for the delta function to give a contribution. But then @(Ep —E^{p)) = 0. 
Altogether the E integration runs from to Ep only: 



3 F 3 F 

N = 2 {^) I dE I ^ 6{E ~ E ' (P)) = 2 (It") / ^ (C - 3) 



=:S(J5) 
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Analogously we calculate the total fermionic energy: 

3 00 

J d 3 p E.(p) = 2 (J^j j d E J d 3 p 5{E - E_(p))e(E- - E.{p))E^{p) 



Kot = 2 ( 



\2tt / 

Fermi 
volume 



Ej? S p 
3 r r /r\3 F 



2 (J^J j dE j d 3 p 5{E - £_ {p))E. (p) = 2 j ' dE E j d 3 p 5(E - E_ (p)) 



V2vr 





R3 R3 



J dEES(E). (C.4) 



Let us evaluate the integral S(E) which can be regarded as the surface of constant energy 
in momentum space. The integration over p\ and P2 may be easily performed by introducing 
polar coordinates r := p\ +p\ and ip: 



S(E) = / *, he - = J d * pS ( E - mK -^--£- + ^JrMAn 

m>a V s ' v '/ 



^ =:-/(P3) 

00 27TOO / „ x OOOO 

U 



y <ip 3 y c/<y3 y dr r 5 ( A - - /(p 3 ) j = 7r y dp 3 J du si A 

— r^ n n \ / — n \ 



2 m ^ 



f(P3) ■ 



-00 N ' -00 

In the last step we changed the integration variable to u = r 2 . We next factor out l/2rriN in 
the argument of the delta function and use the property 

5{ax) = |— -Mx) 



00 00 



to get 

S(E) = 2wm N j dp 3 Jdudl 2m N (A - f(p 3 )) -u\. 

-00 ^ ' 
As the integration over u runs from to 00 only, we have 





and thus 



00 

S(E) = 2vm N J dp 3 e(A-f(p 3 )). 



In order to evaluate S(E) we therefore basically have to find the roots of the function 



,2 J2 



9fe ) := A - f in) - B - m „ - ^ - A + ^^V. (C.5) 
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Since g{p 3 ) is even in p 3 the integration may run from to oo only, bringing a factor 2; 

oo 

S(E) = 4Trm N J d P3 e(g(p 3 )). (C.6) 
o 

Let us analyze the behavior of g{p 3 ). Setting cfipz) = we find its extrema to be located at 

p 3 = , p 3 = ^a 2 -4g 2 AV 2 (C.7) 

(considering only positive values of ^3). We now obviously have to distinguish two cases: 

1. a < 2gAT]: The square root gets imaginary and we only have one extremum at p 3 = 0. 
Since linipg^oo g(p 3 ) = —00 it is a maximum. The function g(p 3 ) hence looks like in fig. 
IC.lal and has, depending on the value of E which shifts the graph vertically, one root 
62 or none. 

g(P3) g(ps) 




P3 



(a) a < 2gAV (b) a > 2gAV 



Figure C.l: Different behavior of g{p 3 ) depending on the value of a 



2. a > 2gAi] : F° r P3 ^ g(ps) now has two extrema; the one at p 3 = is a minimum and 
the other one a maximum. The function thus has the form indicated in fig. IC.lbl and 
has, depending on the vertical shift due to E two (61 and 62)) one (62) or no root. 

We now first deal with case 1. In order to evaluate S(E) of (|C.6j) we have to know from 
which value E\ of E on g(ps) has a root. Below this value g{p 3 ) is never positive and hence 
S{E) = 0. Consequently E\ is the lower integration limit in the integrals of (|C.3|) and (|C.4p 
respectively. By setting g(0) =0 Ei is found to be 

Ei - mN 1 " 2 agAr] (C.8) 

For E > Ei the root 62 of g(p 3 ) lies at 



a 2 



b 2 = b 2 {E) = W ^ + 2m N (E - m N ) + a^2m N (E - m N ) + g\ V 2 , (C.9) 
and we thus find 

SW(E) = Q(E - Ei)4irm N b 2 (E), (CIO) 
the superscript denoting case 1. With (|C.3p we are now able to calculate N: 

3 E ~ F 

i) j dEb2 ^ 



... 3 

N W = 8irm N 
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or 



n 



(l) _ m N 



7T Z 



dE b 2 {E). 



(C.ll) 



The integral may be calculated analytically by changing the integration variable to x = 
2rri]y(E — mjy) + g\rj 2 , leading to an integral of the type f dx ^ x + c\ + c 2 ^fx. Staying on 
the same path as in the previous calculations we would now try to solve this equation for 
Ep and plug it into the expression for E^ ot (which still needs to be evaluated) in order to 
express the latter in terms of n_. However, it does not seem to be possible to do so, since the 
result of the integration contains Ep in a nontrivial way. We therefore will have to do this 
replacement numerically, i.e. for given values of n_ and a we numerically solve (|C.lip for 
Ep (which is possible) and can thus calculate E^ ot . Let us find the corresponding expression 
for -E'tot! picking up f|C4[) : 



3 E r 

E$t = 8nm N U£\ J dE E b 2 (E). 



This integral may also be calculated analytically even though resulting in an even more 
complicated expression. The total fermionic energy density due to n_ finally reads 



-tot 



m N 



7T- 



dE Eb 2 (E). 



(C.12) 



Next, we have to study case 2, i.e. a > Ig^r). For this purpose we have to determine for 
which ranges of E the function g(p$) possesses no, two and one root respectively. The value 
E[ below which g{ps) has no root is found by solving g(p3, m ax) = for E, where P3, m ax is the 
location of the maximum, given by the second equation of (|C.7|) . The result is 

^=™-fg. (C13) 

The next crucial value of E is again where g(0) = 0. This value is of course in turn given by 
Ei of (|C.8p . Hence, for E[ < E < E[ the function g(p$) possesses two roots b% and b 2 , where 



b 1 = b 1 (E) = ] J^ + 2m N (E-m N )-a^2m N (E-m N )+g 2 A r ] 2 , (C.14) 

and b 2 is given by (|C.9p . For E > E\ b 2 is again the only root to deal with. Altogether we 
obtain 



5 (2) (S) = Airm N 



0, E < E[ 



b 2 (E)-h(E), E[<E<Ei. (C.15) 
[b 2 (E), E>El 
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Consequently we have to split the integration in (|C.3p and (|C.4p into two parts: 



F 

dE [b 2 (E) - h(E)] + J dE b 2 {E) 



L 3 



m N 



dE b 2 (E) - J dE h(E) 



( EZ 



n 



(2) _ m N 



7T" 



dE b 2 (E)- J dE bi(E) 



(C.16) 



,(2)- 



8irm]\[ 



L ■ rriN 



It' 



-T 

2vr J 



( El 



F 

dE E [b 2 (E) - bi(E)] + J dE Eb 2 (E) 

Ei 



E, 



J dEEb 2 {E)- J dEEb^E) 



k e I 



( EZ 



(2)- _ m N 
£ tot - n2 



J dE Eb 2 (E) - J dE Ebi(E) I . 



(C.17) 



Now we have to repeat this calculation with E + = E + (p). The only thing that changes 
is the function g(ps) of (|C.5p . which in the case of E + now reads (let us denote it with a 
superscript +): 



g + (j>3) = E-m N 



ri 



Pl+9 2 A V 2 - 



(C.18) 



8m n 2m at 2m, at 

This function always has a maximum at = only, i.e. we simply need to repeat the above 
first case. E\ of (|C.8p will be replaced by 



rp+ , a , a 9Ar] 
= m N + h - . 

8mjy IrriN 



while the root b 2 in (|C.9p becomes 



b%(E) = \/^- + 2m N (E - m N ) - aJ2m N (E - m N ) + g\r] 2 



(C.19) 



(C.20) 



Equations (|C.lip and ()C.12p may again be used with the replacements E F — > Ep~, E\ — > Ef 
and b 2 (E) ^b^(E): 



m N 
n + = — 



+ _ rn N 
£ tot- ^ 2 



-'i 



dE &+(£) 



dE Eb+(E). 



(C.21) 



(C.22) 
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Neutron 



Stars 



During its lifetime a star produces its energy mainly through the fusion reactions 

X H + X H -»• 2 H +e + + 1/ + 1.44 MeV 
2 H + X H -»• 3 He +7 + 5.49 MeV 
3 He + 3 He^ 4 He + X H + X H +12.85 MeV. 



As soon as all the hydrogen is burned, gravitation compresses the helium until other fusion 
reaction start. With increasing diversity of possible nuclear reactions the synthesis of heavier 
nuclei is pushed further until most of the matter consists of Fe, Si and close-by elements. This 
is then the point where the fusion reactions stop and the thermal pressure cannot compensate 
the gravitational pressure anymore. If the star has sufficient mass, the thereupon incipient 
contraction crushes the atomic structure resulting in a mixture of electrons and nuclei, where 
the electrons can be treated as an ideal Fermi gas (i.e. neglecting interactions between them) 
in a first approximation. It is then the Fermi pressure of the electrons that prevents the 
object (called white dwarf) from a further contraction. Let us calculate this pressure by 
taking the electrons as degenerate Fermi gas at T = 0. 

Since we consider free electrons the dispersion relation of one electron is given by (let us 
introduce factors of c and h in this section) 



where m e is the electron mass and p = \p\ the electron momentum. (|2.37p gives the connection 
between the Fermi momentum pp and the electron number density n e : 



At T = all states up to pp are filled and the total energy of the electrons is given by (cf. 



E(p) = \/ m^c 4 + p 2 c 2 , 
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where we have introduced the dimensionless momenta 
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and the function 



X £ ( I x 3 (-t _l 3_~2 



F (1 + ^x 2 F + ...) , 



The pressure P e then follows from the thermodynamical relation dE = TdS — PdV with 
T = 0, leading to P e = —dE e /dV. When taking the derivative of (|D.3P with respect to V we 
have to keep in mind, that xf also depends on V; due to (|D.2p we have = CV -1 / 3 with 
some constant C. So we get 



d . df 8xf q / 7" C it / r 

-/W = 777-lTff = -4V1 + 4777777* = "^V 1 + ^ 



av^ " &c F ^3^4/3 3 y 

leading to 

<9P P mic 5 / x 3 - 



Here we see why only the Fermi pressure of the electrons is important. If the particles are 
non-relativistic (xp <C 1) the Fermi pressure goes like 1/m, i.e. only the light electrons have 
an important contribution. As the electrons get relativistic the star becomes unstable, see 
below. 

We want to express P e as a function of the (for simplicity constant) matter density 

p = an e m N , (D.7) 

where a is the mean number of nucleons (of mass mjv) per electron and where the contribution 
of electrons to the matter density is neglected. Solving (|D.7p for n e and plugging that into 
xf yields 

h ( 3tt 2 P \ 1/3 _ . 

x F = • (D.8) 

- \am N J 



m e c 



Solving this for p and setting xf = 1 we get some characteristic density p c (which will turn 
out to be the typical density of a white dwarf star): 

Since xf oc p 1 ^ 3 the case xp -C 1 in ()D.5p corresponds to p <C p C) while x F 3> 1 means 
p 3> p c - Expanding (|D.6p together with (|D.5p in powers of xf and using ()D.8p we are thus 
led to 

x F « 1 = f^ lP 5 / 3 , p« Pc 

^ 3 U4, ^»i U 2 p 4/3 , p»^ c ' 

with i<Ci and containing the prefactors. In both cases we hence have a so called polytropic 
equation of state, i.e. a relation P = K 1 . A necessary condition for the stability of a star is 
however 7 > 4/3, which follows from an analysis of the total energy of the star when varying 
its radius. Through the so called Lane-Emden function the case p ^> p c (which means that 
the electrons become relativistic) may be connected to the mass of the star as follows: 

M c = ~^M sun , (D.ll) 

where Mc is known as the Chandrasekhar mass limit. Consequently, as M — > Mc the 
star becomes unstable. For a white dwarf consisting of iron, a = 56/26 ~ 2.15 and hence 
Mc = 1.27M sun . Already before reaching this mass limit the reaction 

p + e~^n + v e (D.12) 
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starts to take place (as soon as the Fermi energy of the electrons exceeds the energy that is 
released in the /3 decay, which is about 1.5m e c 2 ), which implies a breakdown of the Fermi 
pressure and a further collapse of the star. A typical white dwarf therefore has the approx- 
imately the mass of our sun and p c of (|D.9|) is indeed the characteristic density of a white 
dwarf. Here we talk about approximately two tons per cubic centimeter, while its radius 
corresponds roughly to that of the earth. 

For a star of mass M > Mq the Fermi pressure of the electrons is not able to prevent the 
matter from further contraction. The denser the material the more protons and electrons are 
converted into neutrons via the reaction (|D.12|) . The neutrinos escape and what is left is an 
object that predominantly consists of neutrons; a so called neutron star. It is then the Fermi 
pressure of the neutrons that keep the object stable. Hence, the above formulas may also be 
used for a neutron star if we replace m e by mjv and set a = 1. The characteristic density of 
a neutron star is therefore 

(D - 13> 

which is about 6 • 10 18 kg/m 3 . Accordingly, the particle number density of the neutrons is 
n = p c /m,N ~ 4 • 10 45 m~ 3 = 4fm~ 3 and is thus of the same order as in an atomic nucleus 
(cf. (|2.57p ). The calculation of the mass and the radius of the neutron star with the methods 
used in the case of white dwarfs is however only possible for densities well below p c . In 
this region the radius of such a dense object can be estimated to be approximately 10 km. In 
order to find out the mass limit one has to go to higher densites, where relativistic corrections 
and other processes become important. Different model calculations led to mass limits for 
neutron stars between 1.5M sun and 3M sun . If the mass of a star reaches this mass limit the 
neutrons become relativistic and the object unstable. Then, the gravitational collapse (most 
likely) cannot be stopped by anything and the star will form a black hole. 
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